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FOREWORD 


This  rsport  dooumants  the  saoond  phase  of  c he  effort  performed  for  the  Nhval 
Sea  Systems  Cfflpmapd  Code  63R  under  contract  N00024-77-C-82&1,  under  Project 
Serial  Nbmber  808H180672,  and  oovers  the  period  from  October  1978  to  October 
1979.  &  includes  all  material  found  in  the  three  Quarterly  Reports,  work  performed 
In  the  fourth  quarter,  connective  material,  the  conclusions  of  the  second  phase 
and  the  approach  for  the  third  phase  of  the  study.  This  study  has  as  lie  goal  the 
examination  of  the  application  of  adaptive  filters  to  passive  eonar  bearing  traddiv. 
Major  steps  toward  this  goal  have  been  realised  via  a  combination  of  statistical 
modeling  and  analysis,  with  verification  by  computer  simulation. 
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1.0  INTRODUCTION 


►  Thl*  la  the  final  report  on  phase  2  of  a  study  oootract  to  determine  the  useful¬ 
ness  of  the  LMS  adaptive  alter  In  a  h oaring  tracker  structure  for  passive  sonar. 
Three  potential  advantages  of  such  a  a  true  hire  over  conventional  trackers  motivated 
the  study,  as  folloars  ^ 

ai  While  ocnvantlooal  trackers  require  a  priori  knowledge  of  the  Input 
power  specta  to  optimise  performance,  the  adaptive  tracker  is 
estimating  the  input  spectra.  This  suggests  that  the  adaptive  tracker 
may  be  insensitive  to  Incorrect  a  priori  assumptions  and  variations 
of  the  input  spectra. 


Stnoe  all  the  correlation  information  between  the  split  array  outputs  is 
contained  in  the  adaptive  filter  weights,  the  potential  eatfcps  to  perform 
broadband  and  narrowband  tracking  simultaneously,  J 


o.  Conventional  trackers  are  essentially  designed  for  stationary  targets, 
modified  for  use  with  dynamic  targets  with  the  addition  of  a  tracking 
loop.  However,  the  adaptive  filter  weights  are  adjusted  iteratively 
(dynamically)  to  satisfy  a  mean  square  error  criteria, 

which  may  provide  advantages  against  moving  targets.  - 

The  fl ret  phase  of  this  study  addressed  the  feasibility  of  such  a  tracker 
using  the  Least  Mean  Square  (LMS)  adaptive  filter  algorithm,  and  developed 
the  tracker  structure  of  Figure  1-1.  hi  this  structure,  the  two  beamfoimed 
spilt  array  outputs  are  provided  as  the  two  inputs  to  an  LMS  adaptive  Alter 
configured  as  a  canceller.  The  output  of  one  half  array  processed  by  a  non- 
recu reive  adaptive  filter,  is  subtracted  from  (he  other  half  array  to  yield  an  error 
signal.  The  error  signal  is  then  used  to  recursively  update  the  adaptive  filter 
weights,  or  impulse  response  such  that  the  mean  square  error  ie  minimised. 

If  x(nTs)  and  d(nT«)  are  the  two  inputs,  with  Ta  the  sample  rate  and  the  time 
Index,  the  adaptive  Alter  stores  the  data  vector 

Xpi)  -  [x£hTa]  x((n-l)Ts]  ....  x£pi-M+l)(r*rT  (1) 

where  XT  denotes  X  transpose.  The  output  of  the  adaptive  filter  is 

y(n)  -  2TpifiC(n)  (1) 

where  W(n)  is  a  vector  of  M  weights  at  time  nTa.  The  error  output  is  therefore 

(pi)  -  dpt)  -  ffTpi)X(n).  (3) 

The  LMS  algorithm  updates  the  weight  vector  on  each  iteration  aa 

HPi«l)  •  W(n)*  m  ( (n)X (n)  (4) 

where  p  la  a  weight  t^data  coefficient.  This  parameter  controls  the  rate  of 
convergence,  algorithm  noise,  and,  ultimately,  the  stability  of  the  algorithm. 


Figure  1-1.  Adepdr*  1 
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k>  converging  to  minimize  the  mean  square  error  between  half  array  outputs, 
the  adaptive  filter  must  incorporate  any  time  delay  (or  for  narrowband  souroea, 
phase  shift)  between  its  inputs  In  the  weight  vector.  For  broadband  Inputs,  the 
weights  have  the  same  shape  as  the  signal  autocorrelation  function  with  the  peak 
located  at  the  delay  between  the  phase  centers.  Hence,  the  tracker  must  determine 
(he  location  of  the  peak  of  the  weights,  using  Interpolation  between  discrete  taps 
If  neosssary.  For  sinusoidal  inputs,  the  weights  oocverge  to  a  sinusoid  with  phsse 
equal  to  the  phase  shift  between  array  halves.  The  tracker  than  determines  the 
phase  shift  by  means  of  spectrum  analysis  of  the  weight  vector.  Either  the  time 
delay  or  phase  shift  la  easily  converted  to  a  bearing  estimate. 

The  first  phase  of  the  study  also  provided  many  of  the  analytical  techniques 
required  for  analysis  of  the  adaptive  tracker,  ta  particular.  If  provided  the 
frequency  domain  model  for  the  adaptive  tracker  abeam  in  Figure  1-2.  b  the 
frequency  domain  version,  the  two  Inputs  are  Fast  Fourier  Transformed  (FFT'd) 
and  a  single  tap,  oomplex  LMS  adaptive  filter  performed  Independently  b  each 
FFT  bln.  It  was  shown  that  this  configuration  performs  equivalently  to  the  time 
domain  filter,  both  In  transient  and  steady  state  operation,  under  certain  condi¬ 
tions  on  the  filter  design  parameters.  The  primary  advantage  of  the  frequency 
domain  model  la  that  both  the  mean  and  variance  of  the  weights  can  be  predicted. 

By  careful  selection  of  filter  parameters,  this  model  can  then  be  used  to  predict 
the  variance  of  the  time  domain  weights,  which  cannot  be  obtained  by  analysis 
in  the  time  domain. 

Given  the  statistics  of  the  time  and  frequency  domain  weights,  it  eras  possible 
to  obtain  predictions  of  the  bearing  estimate  starts  tics  for  broadband  Inputs.  A 
continuous  adaptive  tracker  eras  shown  to  perform  within  0. 5  dB  of  the  Cramer- 
Rao  Lower  Bound  on  the  variance  of  all  unbiased  estimators.  Performance  results 
were  also  derived  for  a  practical  realization  of  the  tracker,  operating  at  a  reason¬ 
able  sampling  rate  with  Interpolation  between  samples  providing  fine  delay  (hence 
bearing)  resolution.  These  results  were  shown  to  be  comparable  to  a  conventional 
tracker.  The  first  phase  of  the  study  also  developed  the  mean  weight  of  the 
adaptive  filter  with  dynamic  broadband  and  narrowband  Inputs.  The  results  are 
applied  here  in  the  analysis  of  tracker  performance  for  dynamic  targets.  Finally, 
the  first  phase  developed  extensive  simulation  programs  for  the  tracker  structure, 
which  were  used  to  validate  the  analyaee,  and  provide  an  initial  look  at  dynamic 
tracking  performance. 

The  first  final  report  established  the  feasibility  of  the  adaptive  traaker  structure, 
particularly  for  broadband  targets.  This  report  addresses  In  greater  detail  the 
potential  advantages  given  in  Items  a  through  c,  above.  The  material  given 
oonalsts  of  the  results  reported  In  the  first  three  quarterly  reports,  plus  work 
done  during  the  fourth  quarter.  Section  2.0  summarises  the  results  of  the  study 
without  detailed  description  of  mathematical  models  or  derivations,  which  are 
included  as  appendices.  Same  expressions  and  figures  from  the  appendices  are 
repeated  In  the  summary  for  convenience.  Section  S.  0  draws  general  conclusions 
and  relates  the  result  given  in  this  report  to  efforts  planned  for  phase  3. 
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2.  0  SUMMARY  OF  RESULTS 


ThU  section  presents  the  results  of  the  seoond  phase  of  the  Adaptive  Tracking 
System  Study  In  summary  form,  without  detailed  derivations  and  analyses,  which 
are  given  In  the  appendices.  The  results  fall  Into  six  categories,  discussed 
separately  below. 

Section  2. 1  oooslders  the  bearing  estimation  performance  of  the  narrowband 
(NB)  adaptive  tracker  with  a  narrowband  target  In  unoorrelated  background  noise 
The  statistics  of  the  bearing  estimate  are  derived,  and  the  variance  compared  with 
the  Cramer-Rao  Lower  Bound.  In  section  2. 2,  the  effects  on  the  NB  adaptive 
tracker  of  broadband  noise  emanating  from  the  same  target  are  studied.  Using 
essentially  the  same  analysis.  Section  2.3  considers  the  NB  adaptive  tracker  with 
a  NB  target  in  the  presence  of  a  broadband  Interference  from  another  direction.  In 
particular,  the  bias  due  to  the  Interference  and  the  effects  on  the  bearing  estimate 
variance  are  given. 

Ctoe  motivation  for  studying  the  adaptive  tracker  has  been  its  possible  Insensi¬ 
tivity  to  Incorrect  assumptions  on  input  spectra  relation  to  conventional  trackers. 
Section  2.4  discusses  the  relative  performance  of  the  adaptive  tracker  and  a 
conventional  split  beam  tracker  (SBT)  when  either  the  noise  or  signal  spectra 
differ  from  those  assumed. 

Section  2. 5  deals  with  the  performance  of  the  broadband  adaptive  tracker  with 
a  broadband  target  and  a  broadband  interference.  Both  a  continuous  realisation  and 
a  discrete  version  of  the  adaptive  filter  are  considered.  Finally,  section  2.6 
considers  the  performance  of  the  broadband  adaptive  tracker  with  a  broadband 
target  whose  bearing  Is  changing  linearly  with  time.  This  Is  analytically  compared 
with  a  Bearing  Deviation  Indicator  (BDI)  with  the  same  dynamic  Input. 

*.  1  Narrowband  Adaptive  Tracker  Performance  With  Stationary  Narrowband 
Target.  The  adaptive  tracker  structure  for  estimating  'Bearing  from  broadband 
and  narrowband  signals  Is  shown  in  Figure  1-1.  The  outputs  of  half  array  beam- 
formers  steered  In  the  vicinity  of  the  target  are  provided  as  inputs  to  aa  LMS 
adaptive  filter  oonflgured  as  a  canceller.  The  Alter,  in  minimising  the  mean 
square  error  between  the  filter  output  and  the  reference  inputs,  would  Ideally 
Insert  a  time  delay  equal  to  the  time  delay  between  half  array  phase  centers. 

This  can  be  Interpreted  as  a  phase  shift  for  narrowband  signals.  This  Ideal  result 
is  corrupted  by  background  noise,  finite  observation  time,  the  discrete  (In  time) 
nature  of  the  filter,  and  correlation  properties  of  the  signal  and  noise. 

For  a  narrowband  target,  the  mean  adaptive  filter  weights  converge  to  a 
sinusoid  at  the  frequency  of  the  target  line,  with  a  phase  equal  to  the  phase  shift 
between  array  halves.  To  extract  this  phase,  the  weight  vector  is  Fast  Fourier 
Transformed  (FFT'd),  and  the  phase  of  the  FFT  bln  containing  the  signal  is  used 
as  the  phase  estimate.  Conversion  of  the  phase  estimate  to  bearliy  requires 
knowledge  of  the  line  frequency,  which  is  estimated  as  the  center  frequency  of  the 
FFT  bin  containing  the  signal,  assumed  to  be  the  one  with  the  largest  magnitude. 


r 


During  Phase  1  of  this  study  a  frequency  domain  model  for  the  adaptive 
tracker  was  developed  which  allowed  determination  of  both  the  mean  and  varlanoe 
of  the  adaptive  tracker  weights  for  broadband  Inputs-  Using  model,  the  mean 
and  variance  of  the  FFT  of  the  weight  vector  for  the  case  of  a  narrowband  target 
in  broadband  noise  are  derived,  as  given  In  the  first  section  of  Appendix  A.  Based 
upon  the  large  number  of  Iterations  required  for  convergence  of  the  adaptive  filter, 
the  frequency  domain  weights  are  assumed  to  be  complex  gausslan  random  variables. 
This  allows  determination  of  the  density  of  the  phase  in  each  FFT  bln  by  transforma¬ 
tion  to  polar  coordinates,  then  integrally  out  the  magnitude.  The  bearing  is 
related  to  the  phase  by 

®  (8) 

where  c  is  the  speed  of  sound,  d  the  distance  between  half  array  phase  centers, 
e  the  phase  estimate  In  the  FFT  bin,  and  u  the  estimated  radian  frequency  of  the 
line,  approximated  here  by  the  FIT  bln  oentSr.  By  numerical  integration  of  this 
function  and  Its  square  over  the  phase  density,  the  mean  and  variance  of  the  bearing 
estimate  are  obtained.  The  analysis  assumes  that  the  bln  containing  the  is 
selected  for  phase  extraction.  This  means  that  either  the  signal -to- noise  ratio  In 
the  bln  must  be  sufficiently  high  or  else  some  other  means,  such  as  aAiiH«n«i  non¬ 
coherent  integration  must  be  used  to  Identify  the  bln  containing  the  aigml  if  the 
predictions  are  to  be  valid.  The  details  of  the  analysis  are  given  In  Appendix  A. 

The  analysis  shows  that  the  bearing  estimate  is  essentially  unbiased,  and  this 
Is  confirmed  by  simulations.  Figure  2-1  shows  the  r.  m.  s  bearing  error,  determined 
as  above,  for  a  particular  adaptive  tracker.  Tbs  signal  Is  a  100  Hx  sinusoid 
emanating  from  a  target  at  0°  bearing  and  observed  by  an  array  with  150  feet 
between  half  array  phase  centers.  The  FFT  Is  an  128  point  transform  with  2  Hi 
frequency  resolution.  Also  shown  In  the  figure  are  the  results  of  computer  simula¬ 
tions  of  the  same  tracker,  demonstrating  excellent  agreement  with  the  analytical 
predictions. 


An  approximation  for  the  variance  of  the  phase  In  the  FFT  bln  containing  the 
signal  la  derived  using  a  Taylor  expansion  of  the  Inverse  tangent,  as  used  to 
obtain  the  phase  from  the  FFT  weights.  By  retaining  only  second  order  lower 
terms,  then  Integrating  the  result  over  the  phase  density,  It  Is  possible  to  obtato 
simple  approximations  to  the  phase  variance  at  high  signal  to  noise  ratio.  For 
bearing  angles  near  the  maximum  response  axis  (MBA)  of  the  array,  this  Is  easily 
mapped  to  s  bearing  variance  approximation,  given  by 


Tar  [8]  - 


M2^ 

(2*k)2 


3>  4>  ♦  4Y  ♦  1 


,2  ' 
m 


>  2(V«-1)  - 


(y2*4y*: 


-] 


(6) 
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where  M  is  the  number  of  FFT  points,  k  the  Index  of  the  FFT  bin  containing  the 
signal.  T,  the  sample  rate,  >  the  signal  to  noise  ratio  In  the  FFT  bln,  andrT 
the  adaptive  filter  time  constant  in  iterations.  This  can  be  further  approximated 


6 


Figure  2-1:  Adaptive  Tracker  RMS  Bearing  Error  For  Static 
Narrowband  Target. 
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Hence,  at  high  SNR,  the  variance  1 *  Inversely  proportional  to  the  SNR  and  the 
“Integration  time”,  r.,  within  the  filter.  The  approximation  given  In  (6)  la  ahown 
plotted  with  the  exact  numerical  evaluation  for  an  adaptive  tracker  In  Figure  2-2. 
The  tracker  configuration  la  the  tame  aa  la  Figure  2-1.  Over  the  range  In  which 
the  numerical  analyai*  la  valid,  the  approximation  la  excellent. 


In  Appendix  B,  the  Oamer-Rao  Lower  Bound  (CR  LB)  for  the  narrowband 
tacking  problem  la  derived.  The  CR  LB  la  a  lower  bound  on  the  variance  of  any 
unbiased  estimate  given  the  same  set  of  observations.  The  observations  are  taken 
to  be  N  FFTs  of  the  split  array  beam  former  outputs.  By  using  the  approximation, 

(7),  from  above,  it  la  shown  that  at  high  SNR,  the  adaptive  tracker  bearing  variance, 
is  0.  5  dfi  above  the  bound.  Plots  for  the  adaptive  tracker  configuration  given  above 
for  the  plot  of  Figure  2-1  Indicate  that  for  rma  bearing  error  lees  than  1°,  the 
adaptive  tracker  requires  1  to  5  dB  more  SNR  to  achieve  variance  equal  to  the 
CRLB.  It  should  be  noted  that  the  CRLB  la  not  necessarily  attainable  by  any 
estimator.  Figure  2-3  la  typical  of  the  performance  of  the  adaptive  tracker  relative 
to  the  CRLB. 


2.2  Narrowband  Adaptive  Tracker  Performance  Wltfa  A  Combined  Broadband 
And  Narrowband  Target.  Tha~sams  general  approach  used  In  the  preceding  section 
la  also  employed  In  determining  the  performance  of  the  narrowband  adaptive 
tracker  with  a  target  rvdtadag  both  narrowband  and  broadband  components.  The 
only  difference  In  the  analysis  la  the  determination  of  the  mean  and  the  variance 
of  the  adaptive  tracker  weights  with  this  Input.  In  Appendix  C  the  mean  and  var¬ 
iance  of  the  arelght  vector  la  rederived  with  a  narrowband  directional  source  and 
a  broadband  directional  source,  not  necessarily  (Tom  the  same  direction.  When 
the  two  directions  are  set  equal,  the  stall  sties  are  those  required  here.  The 
broadband  and  narrowband  components  are  treated  aa  independent  in  the  analysis. 
The  broadband  signal  gives  rise  to  cross  terms  between  the  broadband  and  nar¬ 
rowband  components  in  the  variance  expression,  aa  well  aa  broadband  compo¬ 
nents  In  both  the  mean  and  variance. 

Procedi ng  as  In  Section  2. 1,  the  weight  statistics  are  mapped  to  bearing 
statistics  under  the  assumption  that  the  weights  are  gaussian,  by  means  of 
coordinate  transformation  and  numerical  Integration  ae  detailed  in  Appendix  D. 
ft  Is  shown  thst  when  the  narrowband  component  Is  not  centered  In  the  FFT  bin, 
a  bias  In  the  bearing  ee tlm ate  results.  From  Appendix  D,  this  bias  Is  approximately 
given  by 


(«) 
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(civil* 


when 


-  actual  target  bearing 
wQ  •  narrowband  line  frequency 
k  •  Index  of  FFT  bln  containing  algnal 
M  -  FFT  alas 
Tg  -  aample  rate 

Thla  blaa  decreases  aa  the  FFT  reaolutlon  ta  made  finer.  Since  the  max¬ 
imum  error,  given  adequate  SNR  to  select  the  bin  nearest  the  signal,  la  one- 
half  of  an  FFT  bln,  the  bias  la  proportional  to  the  ratio  of  FFT  reaolutlon  to 
center  frequency. 


where 

) _ -  largest  bearing  angle  considered 

m  max 

tlR  -  FFT  resolution 

Thla  relationship  can  be  used  to  set  the  FFT  resolution  for  acceptable  bias. 
Note  that  the  estimate  Is  unbiased  for  targets  at  broadside. 

The  rms  bearing  error  of  a  narrowband  adaptive  tracker  with  a  composite 
target  la  shown  In  Figure  2-4  as  a  function  of  the  narrowband  algnal  to  background 
noise  ratio  for  various  levels  of  broadbsnd  signal.  The  presense  of  broadband  noise 
from  the  same  target  can  significantly  improve  the  narrowband  estimate,  particu¬ 
larly  when  its  power  exceeds  that  of  the  narrowband  component.  The  tracker  here 
Is  a  128  tap  filter  with  FFT  resolution  of  2  Ha  and  a  weight  update  ooeffidma,  u , 
of  2~1*.  The  array  has  160  ft  between  phase  centers,  and  the  target  Is  at  0° 
bearliv  with  a  100  Ha  Une  and  a  flat  broadband  spectrum.  As  In  section  2. 1,  these 
results  are  valid  only  when  the  correct  FFT  bln  la  selected.  Also  shown  In  the 
figure  are  the  results  of  computer  simulations,  which  verify  the  analytical 
predictions. 


By  noting  the  performance  of  the  traoker  In  Figure  2-4  at  low  narrowband 
SNR  for  moderate  to  high  broadband  SNR,  that  la  whan  moat  target  energy  la 
broadband.  It  appears  that  the  narrowband  traoker  la  a  good  estimator  of  the 
bearing  of  a  broadband  target.  However,  the  analyses  here  assume  that  the  delay 
between  phase  centers  is  always  an  Integer  number  of  aample  intervals,  aa  assump¬ 
tion  necessary  In  the  uncorrelated  notes  anas.  For  oorrelated  noise  and  non-integer 
delays,  this  may  change  somewhat.  Further  analysis  and  simulation  would  be 
required  to  conclude  that  the  narrowband  tracker  can  always  be  used  for  broad¬ 
band  signals. 
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Figure  2-4:  RMS  BMrtaf  Error  For  Narrow  bond  Adaptive  Tracker  And  Target 
With  Broadband  And  Narrowband  Radiated  Nolee. 


♦C  <tft« 


The  approximation  developed  to  Appendix  A  for  the  nna  bearing  error  wltb2 
a  narrowband  signal  only,  (A -46),  la  equally  valid  here  when  the  parameter,  fi  , 
la  replaoed  by  (D-6)  from  Appendix  O,  although  die  complexity  of  the  latter  makes 
the  approximation  lees  useful. 

2.3  Performance  Of  Tim  Narrowband  Adaptive  Tracker  With  A  Narrowband 
Target  And  firoadband  Ihterferenoe.  The  analyses  and  almulationa  In  section  1.1 
and  2.2  consider  a  target  In  a  noise  Held  that  Is  uncorrelated  between  array 
halves.  A  realistic  environment  oould  be  expected  to  Include  directional 
other  than  the  target,  referred  to  as  Interferences.  Appendix  D,  using  d 
results  of  Appendix  C,  considers  the  effects  of  a  single  broadband  lnterf 
the  performance  of  a  narrowband  tracker  with  a  narrowband  target.  The 
statistics  developed  In  Appendix  C  represent  this  situation  when  the  bear 
of  the  broadband  and  narrowband  components  are  not  equal.  The  mapping 
weight  statistics  to  bearing  Is  Identical  to  that  described  in  section  3. 2. 
however,  differ  significantly. 


When  a  broadband  Interference  of  moderate  strength  relative  to  the 
target  Is  present,  the  dominant  effect  la  a  bias  In  the  bearing  estimate  U 
interference.  For  angles  within  ±10°  of  broadside,  this  bias  Is  given  In  i 
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ig  -  target  bearing 
f  j  -  Interference  bearing 
X  *  signal  to  Interference  ratio  In  FFT  bln 
«**()  *  target  frequency 
k  *  Index  of  FFT  bln  nearest  the  signal 
d  ■  distance  between  half  array  phase  centers 
c  *  speed  of  sound 
M  •  FFT  also 
Tg  -  sample  rate 

The  first  term  la  due  to  non-bin  centering  of  die  signal,  and  will  typically  be 
negligible  In  oomparlsoo  to  the  seoond  term.  The  bias  can  also  be  computed 
numerically  by  integrating  over  the  phase  density,  as  described  in  section  2. 1,  and 
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in  Appendix  0.  When  thin  la  dona*  curve*  like  those  shown  in  Figure  2-6  result. 
This  figure  shows  the  bias  Is  a  particular  narrowband  adaptive  tracker  as  a  ftinodon 
of  the  slgnal-to-nolae  ratio  (SNR)  of  a  target  at  0°,  with  an  Interference  at  4°  for 
four  Interference- to-nois*  ratio*  (INK).  Whan  the  SNR  la  low  In  comparison  to  the 
INR  In  the  FFT  bln,  the  estimate  la  very  nearly  the  Interference  angle  while  at  high 
SNR,  the  estimate  la  essentially  unbiased.  The  results  of  computer  simulations  are 
also  shown  la  Figure  3-6,  and  show  excellent  agreements. 

This  result  should  not  be  considered  as  unique  to  the  narrowband  adaptive 
tracker,  but  will  be  typical  of  tracker*  of  this  generic  type  that  extract  phase 
information  from  FFT  outputs.  This  suggests  that,  at  least  In  part,  the  disappoint¬ 
ing  performance  experienced  with  many  existing  narrowband  trackers  may  be  due 
to  the  comm ulatlve  bias  effect  of  numerous  broadband  directional  sources  such  as 
own  ship  and  distant  shipping. 

The  variance  of  the  tracker  with  a  broadband  interference  bee  also  been 
determined  using  the  numerical  Integrations  of  the  phase  density  described  in 
Appendix  D,  and  the  results  for  s  particular  example  are  shown  In  Figure  2-6. 

The  signal  Is  s  100  Ha  line  received  by  an  array  with  160  ft  between  phase  centers, 
and  emanating  from  a  target  of  4°  bearing.  The  alible  broadband  interference  is 
located  at  0°.  The  128  tap  adaptive  tracker  uses  an  FFT  resolution  of  2  Ha  and  a 
weight  update  coefficient  of  2“**.  The  standard  deviation  of  the  bearing  error  Is 
plotted  as  s  function  of  SNR  for  four  values  of  Inter fcrence-to-notae  ratio.  Any 
evaluation  of  tracking  performance  must  consider  both  Figure  2-6,  the  bias,  and 
2-6,  the  standard  deviation,  simultaneously.  For  example,  while  a  -*-10  dB  inter¬ 
ference  gives  very  low  standard  deviation  at  low  SNR  from  Figure  2-6,  the  estimate 
1*  severely  hissed,  from  Figure  2-5.  A  reasonable  performance  measure  is  the 
root  mean  square  error,  given  by  the  square  root  of  the  sum  of  the  square  of  the 
bias  and  the  variance.  It  is  Interesting  to  note  that  at  high  SNR,  when  the  bias  la 
negligible,  the  Interference  has  the  effect  of  Increasing  the  variance  of  the  estimate. 
Results  of  computer  simulations  of  this  adaptive  tracker  are  also  shown  on  Figure 
2-6,  and  verify  the  analytical  predictions. 

2.4  Performance  Of  The  Broadband  Adaptive  Tracker  With  faoorrect  A  Priori 
Aaamnptlooa  Cto  Input  Spectra.  The  broadband  adaptive  tracker  has  been  suggested 
ss  an  alternative  to  conventional  split  beam  correlator  tracker*  because  a f  several 
potential  advantages  stemming  from  the  adaptive  characteristics.  One  of  these 
possible  advantages  la  that  the  adaptive  filter  does  not  require  a  priori  knowledge 
of  the  signal  and  noise  statistics,  but  estimates  these  dynamically  during  operation. 
Conventional  split  beam  trackers,  on  the  other  hand,  require  the  Input  signal  and 
noise  spectra  in  order  to  optimise  their  performance  (see,  for  example,  [2]  or  BJ. 
The  goal  of  this  analysis  is  to  determine,  at  least  for  some  simple  oaaee,  whether 
or  not  the  adaptive  tracker  Is  more  tolerant  of  spectral  variation  about  those 
assumed  a  priori  than  the  split  beam  tracker  (SBT). 

The  approach  la  as  follows.  A  signal  and  noise  spectrum  la  assumed  a  priori, 
and  the  weighting  filter  for  the  SBT  is  set  using  these  spectre.  Then,  with  the  Input 
spectra  matching  the  a  priori  assumption*  the  variance  of  the  bearir«  estimate  for 
the  adaptive  tracker  and  SBT  ars  computed  analytically.  Either  the  input  noise  or 
signal  spectrum  la  then  varied  without  ahaeytng  the  weigh  tlx*  Alter  of  the  SBT,  and 
the  variances  recomputed.  The  ratio  of  the  variances  with  correct  and  Incorrect 
s  priori  knowledge  gives  the  change  In  performance  of  the  two  trackers.  These 
are  then  used  to  determine  which  of  the  two  is  less  susceptible  to  the  rhsigflng 
inputs. 
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Figure  2-6:  RMS  Bearing  Error  For  The  Narrowband  Adaptive  Filter  hi  the 
Preeeoee  Of  A  Broadband  Interference. 
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to  Appendix  4  of  Q],  the  mean  and  ▼arlanoe  ol  the  transfer  function  of  a 
continuous  adpetlve  filter  waa  derived  under  the  condition  of  equal  signal  and  nolae 
banda.  Thla  waa  done  by  reaolvfng  the  Input*  Into  deaerate  frequency  bln*,  and 
applying  a  single  tap  oomplex  adaptive  filter  In  each  bin,  aa  dlaonaaed  In  Appendix 
9  of  DL3.  The  tranafer  functlona  are  obtained  aa  the  limit  of  the  dlacreta  oaae.  Thla 
approach  depeoda  upon  the  independence  of  the  Individual  frequency  bina.  The 
reaultlng  tranafer  function  can  be  uaed  to  obtain  the  Impulae  reaponae  of  the  adaptive 
filter.  h(r). 

The  angle,  9  ,  la  to  be  eetlmated  by  determining  the  value  of  r  for  which  the 
Impulae  response,  h(r),  la  a  maximum.  Thla  la  an  eatlmate  of  the  delay  between 
arrtvala  of  the  signal  wavefront  at  the  half  array  phase  centers,  and  can  be  con¬ 
verted  to  a  bearing  estimate,  9  ,  using 

r  -4«ln  9  .  CL1) 

c 

The  value  of  9  for  which  h(r)  la  peaked  oorreepooda  to  the  value  for  which  the 
derivative  of  h(r>  with  respect  to  9  1*  aero.  In  the  neighborhood  of  the  aero  cross¬ 
ing  the  derivative  of  the  Impulse  response  1*  approximately  linear.  The  fluctuations 
in  the  derivative  of  the  transfer  function  then  map  through  the  linear  function  to 
provide  an  estimate  of  the  mean  square  error  In  the  angle  estimate.  Thu*  the  mean 
value  of  the  derivative  of  the  transfer  function,  the  standard  deviation  of  the  deriva¬ 
tive  of  the  tranafer  function,  and  the  slope  of  the  mean  of  the  derivative  of  the 
tranafer  function  are  required,  at  the  point  where  9=9.  Then  the  errors  In  the 
transfer  function  can  be  mapped  Into  the  errors  In  one's  ability  to  extract  the  peak 
of  the  tranafer  function  aa  foUowe : 
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This  approach  ia  motivated  by  the  treatment  of  BDI  split  beam  trackers  by 
MacDonald  and  Schulthelss  In  Reference  [2]. 

to  Appendix  E,  the  adaptive  tracker  la  compared  witothe  conventional  split 
beam  tracker,  »»«"g  and  Schulthelss '  results.™  These  results  require 

a  priori  assumptions  regarding  the  signal  and  nolae  spectra.  The  a  priori  signal 
and  nolae  spectra  are  denoted  S.a(w)  and  S^lu),  respectively,  not  necessarily 
equal  to  the  actual  spectra. 

The  effects  of  incorrect  a  priori  assumptions  on  the  performance  of  the  split 
beam  tracker  are  examined  in  several  specific  cases.  This  la  done  by  comparing 
the  variance  of  toe  bearing  estimate  under  actual  conditions  with  the  oeae  whan  a 
priori  assumptions  arc  true.  The  perform  an  oe  of  the  adaptive  Alter  under  both 
conditions  la  also  computed,  to  order  to  maintain  toe  same  time  constant  for  the 
adaptive  tracker,  the  power  at  the  input  la  normalised  to  unit  power  (equivalent 
to  an  AGC). 
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The  two  special  oases  shown  In  Figure  2-7  are  studied  in  detail.  In  case  1, 
the  signal  and  noise  spectra  are  assumed  a  priori  to  be  Oat  over  a  bound  0  to  u 
with  levels  Le  and  Lq,  respectively.  The  actual  noise  spectrum  is  as  assumed,* 
but  the  actual  signal  spectrum  differs  in  level  from  the  assumed  In  the  region 
from  u,  to  The  analysis  is  performed  at  low  slgnal-to-nolae  ratio.  It  Is 
shown  in  Appeodtx  E  that  at  low  SNR,  the  performance  of  the  SBT  and  the  adapt¬ 
ive  tracker  changes  in  the  same  way  with  changing  signal  spectrum.  This  Is 
due  primarily  to  the  low  SNR,  since  the  adaptive  filter  is  responding  mostly  to 
the  noise  Input  under  the  minimum  mean  square  error  criterion. 

The  second  case  shown  In  Figure  2-7  uses  the  same  a  priori  assumptions  as 
case  1,  but  the  actual  noise  spectrum  differs  from  the  assumed,  and  the  signal 
spectrum  matches  the  assumptions.  The  noise  spectrum  differs  in  level  In  the 
region  Mq  to  wa.  hi  this  oase,  the  analysis  of  Appendix  E  shows  considerable 
difference  between  two  trackers.  As  the  noise  decreases  In  the  region  wn  to  ua, 
the  performance  of  both  trackers  improves,  but  the  adaptive  tracker  Improves 
more  rapidly  since  It  can  adapt  to  better  take  advantage  of  the  reduced  noise.  As 
the  noise  Increases,  both  trackers  degrade,  and  except  for  very  small  changes  in 
level,  the  adaptive  tracker  degrades  more  slowly. 

The  two  trackers  can  be  compared  In  terms  of  a  performance  ratio,  D,  given 

by 


Var  (0)/VQ 
D4(  9  )/Dq 


<13) 


where 


Var( »  )  ■  variance  of  adaptive  tracker  estimate  with  actual  Input 
spectra 

V  ■  variance  of  adaptive  tracker  estimate  If  a  priori  assumptions 
were  true 

_ A 

D  (9  )•  variance  of  SBT  estimate  with  weighting  filter  set  using  a 
priori  assumptions,  but  actual  spectra 


D  -  variance  of  SBT  estimate  with  weighting  function  for 
assumed  spectra,  if  assumptions  were  true 


The  numerator  and  denominator  represent  the  change  In  performance  as  the  noise 
spectrum  varies  from  the  assumed  (note  that  when  either  Is  greater  than  unity,  the 
performance  Improves  with  the  change  in  spectrum).  Therefore,  when  D>1,  the 
split  beam  tracker  has  responded  more  favorably  to  the  change  (1.  e. ,  it  has  de¬ 
graded  less  or  Improved  more),  while  D<1  favors  the  adaptive  tracker.  Figure  2-8 
shows  D  plotted  as  a  function  of 


Kn 


actual  noise  level  in  u  tom 

_  n  a 

assumed  noise  level  in  w  tow 

Q  A 
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Thin  la  plotted  for  several  values  of  K  ,  where 

v 


w 


la  a  measure  of  the  site  of  the  region  in  which  the  assumptions  are  invalid.  As 
men  Honed  above,  for  Kq<  1  and  for  K-  significantly  greater  than  1,  the  adaptive 
tracker  yields  superior  performance  (D<  1).  When  the  noise  level  differs  signifi¬ 
cantly  from  that  assumed,  the  advantage  of  the  adaptive  tracker  can  be  great. 

For  small  Increases  In  the  noise  level,  the  SBT  shows  a  slight  advantage  over 
the  adaptive  tracker,  as  shown  for  Kq  slightly  greater  than  one  in  Figure  2-8.  A 
Is  shown  In  Appendix  C  that  the  advantage  1s  at  most  on  the  order  of  0. 5  dB.  This 
must  be  compared  with  the  much  larger  advantages  of  the  adaptive  filter  for  other 
values  of  Kq.  Further,  it  should  be  noted  that  when  <  o.  5  the  a  priori  assumptions 
as  to  spectra  shape  are  improving  with  decreasing  Kw.  The  worst  case  situation 
for  spectral  shape  Is  then  Kw  -  0. 5.  With  this  In  mind,  the  appendix  shows  that 
increases  In  the  noise  up  to  about  3  dB  favor  the  SBT,  but  beyond  that,  the  adaptive 
tracker  is  superior.  The  adaptive  tracker  appears  then  to  offer  significant  Immunity 
to  changing  input  spectra  relative  to  the  SBT  for  roost  spectral  variations  studied, 
which  the  SBT  offers  only  a  slight  relative  advantage  when  the  change  Is  a  small 
Increase  In  noise. 


2.5  Performance  Of  The  Broadband  Adaptive  Tracker  In  The  Presence  Of 
A  Broadband  Interference.  The  flrst  phase  ofthe  Adaptive  Tracking  System  Study 
considered' the  performance  of  the  broadband  adaptive  tracker  with  a  single  broad¬ 
band  directional  source,  the  target.  Cl]  This  was  done  by  using  the  frequency 
domain  model  of  the  tracker  to  develop  expressions  for  the  mean  and  varlanoe  of  the 
adaptive  filter  weights.  These  frequency  domain  weights  were  then  transformed 
to  give  the  second  order  statistics  of  the  time  domain  weights,  or  Impulse  response 
of  the  filter.  The  delay  between  half  array  phase  centers  was  extracted  from  the 
impulse  response  by  determining  the  location  of  the  maximum.  The  weight  vector 
statistics  must  therefore  be  mapped  into  the  mean  and  variance  In  the  determination 
of  the  peak. 


This  was  done  using  the  method  used  by  McDonald  and  Schulthelsa .  ^  The 
value  of  the  peak  ooi  reeponds  to  the  value  for  which  the  derivative  of  the  impulse 
response  with  respect  to  bearing  Is  aero.  Near  the  aero  crossing,  the  derivative, 
h  ’  (t>,  is  approximately  linear.  The  fluctuations  in  the  derivative  than  map  through 
the  slope  of  the  derivative,  3[h’(t)]  3»,  to  the  fluctuations  In  the  aero  crossing, 
or  peak  selection  process.  The  variance  in  the  bearing  estimate  la  then 


Vtr1  *[•  ] 


(14) 


where  hfr)  is  the  Impulse  response  of  the  sdsptive  tracker  and  t*  Is  the  location 
of  the  aero  crossing. 


Two  cams  wore  considered,  a  continuous  realisation  of  the  adaptive  Alter, 
and  a  discrete,  or  sampled,  vers  loo.  The  continuous  case  was  analysed  by  deter¬ 
mining  the  mean  and  covariance  functions  of  a  continuous  frequency  Hrenein  tracker 
aa  the  Umlttlng  case  of  the  discrete  tracker.  Fourier  transforms  than  ylslded 
the  statistics  of  the  Inpulse  response  required  to  evaluate  (14).  The  rwults  of  the 
analysis  was  primarily  of  value  as  a  means  of  comparison  with  the  Cramer- Rao 
Lower  Bound  (CRLB)  an  the  variance  of  any  estimator  with  the  same  observation. 

In  practice  however,  the  adaptive  tracker  la  discrete  In  time,  and  Interpolation 
must  be  used  to  determine  the  location  of  the  peak  between  discrete  tape.  The  Inter¬ 
polation  significantly  increases  the  variance  above  that  obtained  a 
adaptive  filter.  This,  then.  Is  the  more  realistic  performance  measure  from  a 
practical  point  of  view. 

Appeodtx  F  rap  seta  the  analysis  of  the  variance  of  the  beating  aatlmate  of  a 
broadband  target,  but  includes  a  single  whits,  broadband  Interference  at  another 
bearing.  The  derivation  of  the  mean  and  variance  of  the  sd^>Uve  filter  weights 
usas  the  frequency  domain  model  and  the  results  of  Appendix  V  of  DU  which  gave 
the  mean  and  variance  of  the  frequency  domain  weights.  The  weight  vartanoe 
la  of  the  ume  form  as  In  the  absence  of  Interference,  but  with  the  signal  power 
replaced  by  the  oombtned  signal  plus  Interference  power. 
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where 

Wk*~>  "  •***•  weight  In  k1^  FFT  bln 

M  -  number  of  adaptive  filter  tape 

ft  -  weight  update  coefficient 
2 

<rQ  -  noise  power 

2 

<r  -  signal  power 
o-j  -  interference  power 
The  mean  steady  stats  weight  in  the  k^  bln  la 
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where  A  and  >  ere  the  delays  between  half  array  phase  centers  expressed  In  units 
of  the  sample  Interval  for  signal  and  Interference,  respectively. 


This  expression,  (16),  points  out  the  primary  difference  between  the  oases  of 
interference  and  no  Interference.  When  (16)  is  Inverse  transformed,  it  yellds  the 
sum  of  two  uaimodal  functions,  with  their  peaks  located  at  t  •  AT,  and  t  ■VTg. 

If  the  signal  and  interference  correlation  functions  are  narrow  In  comparison  to  the 
distance  between  the  peaks,  the  resulting  adaptive  tracker  lnpulse  response  will  be 
blmodal,  with  both  peaks  biased  somewhat  toward  each  other.  In  this  case,  if  it 
la  possible  to  identify  the  signal  peak,  either  through  previous  tracking  history, 
or  a  priori  knowledge,  it  will  be  possible  to  obtain  a  relatively  unbiased  estimate 
of  target  bearing,  although  the  variance  will  be  increased  by  the  interference  power. 
Cto  the  other  hand.  If  the  correlation  functions  are  wide  relative  to  the  distance 
between  the  peaks,  the  resulting  weight  vector  may  be  uni  modal,  with  the  peak 
somewhat  between  the  correct  delays  corresponding  to  signal  and  Interference,  hi 
this  case,  only  a  biased  estimate  is  possible. 

expressions  for  the  variance  of  the  adaptive  tracker  bearing  (or  delay)  estimate 
are  developed  In  Appendix  F  as  a  function  of  the  location  at  which  the  tero  crossing 
occurs.  This  is  done  both  for  the  continuous  and  the  discrete  adaptive  filter.  In 
both  cases,  the  expressions  are  somewhat  more  complicated  than  the  target  only 
case,  and  yield  little  insight  without  numerical  evaluation.  Further,  it  la  not 
possible  to  determine  the  location  of  the  peak  or  peaks  In  the  lnpulse  response 
explicitly.  This  must  be  done  numerically  on  the  computer,  then  the  resulting 
location  of  the  peak  substituted  into  the  variance  expressions  for  evaluation.  The 
appendix  presents  the  required  expressions  for  evaluation,  but  numerical  computa¬ 
tion  of  the  results  has  not  been  completed.  Tills  will  be  carried  out  es  part  of 
Phase  3  of  the  study. 


2.6  Adaptive  Tracking  Of  Moving  Broadband  Targets.  As  mentioned  in  the 
Introduction,  one  of  the  reasons  for  studying  the  adaptive' tracker  structure  of 
Figure  1-1  Is  that  it  seemed  possible  that  the  dynamic  adjustment  of  the  tracker 
weights  to  a  minimum  mean  square  error  criterion  would  yield  same  advantage 
over  conventional  tracker  when  the  target  is  moving.  Appendix  G  considers  this 
possibility  by  analytically  comparing  the  adaptive  tracker  and  a  conventional  firet 
order  split  beam  tracker  with  a  broadband  target  at  low  SNR  whose  bearing  is 
changing  linearly  with  time. 


Using  the  expression  for  the  mean  weight  vector  at  low  SNR  with  s  moving 
target  developed  In  Appendix  VI  of  reference  dX  and  assuming  an  exponential 
correlation  function  for  the  tracker  signal  Input,  it  is  possible  to  develop  closed 
form  expressions  for  the  mean  weights  of  the  adaptive  tracker.  These  results 
are  derived  In  Addenda  GI  and  GH  of  Appendix  G.  By  taking  limits  to  obtain  steady 
state  results,  then  determining  the  location  of  the  peak  of  the  mean  weight  vector, 
an  expression  for  the  lag  of  the  adaptive  Is 
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where 

c  -  linear  rate  of  change  of  delay  between  half  array  phase  centers 
o*  -  signal  power 
<r*  -  noise  power 

D 

p  ■  e"b  -  exponential  correlation  function  of  signal 

b  •  bandwidth  of  signal 

6  •  algorithm  sample  time 

n  -  weight  update  coefficient  of  adaptive  filter 

When  »  bed,  then  this  Is  approximately  given  by 
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In  sections  2.4  and  2. 5,  and  In  Appendix  IV  of  [11  a  method  of  determining 
the  variance  of  the  delay  or  bearing,  estimate  from  the  mean  and  variance  of  the 
adaptive  tracker  weights  was  described.  This  mapping,  used  by  McDonald  and 
Scbulthelss  0®]  in  their  analysis  of  conventional  spilt  beam  trackers,  gives  the 
estimate  variances  as 


<1»> 


with  h(t)  the  adaptive  tracker  weight  function  and  tp  the  location  of  the  peak.  To 
obtain  the  statistics  of  the  derivative  required  by  (191,  It  Is  necessary  to  have 
not  only  the  mean  weight  vector,  but  its  variance.  Appendix  G  assumes  small 
slgnal-to-nolse  ratio,  and  that  the  variance  does  not  differ  significantly  Grom 
the  noise-only  case  at  low  SNR.  The  variance  used  in  (19)  is  therefore  the 
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’•null  hypothesis"  variance.  With  thia  approach  the  variance  in  the  delay  esti 
mate  is  found  to  be 


This  result  screes  with  previous  results  P3  for  (he  static  case  when  c  -»0, 
and  demonstrates  the  expected  behavior  with  SNR,  rate  (c),  and  bandwidth. 

The  lag  variance  of  the  conventional  split  beam  tracker  are  determined 
■ii|pg  a  linearized  mathematical  model  described  in  Appendix  G.  Using  the  results 
of  this  analysis,  the  design  parameters  of  the  conventional  tracker  and  the  adaptive 
tracker  are  selected  so  thst  the  two  have  equal  lags  for  a  particular  dynamic  al 
input.  This  allows  the  two  trackers  to  be  compared  in  terms  of  the  delay  variance 
for  equal  lays.  The  ratio  of  the  variance  in  the  delay  estimate  from  the  adaptive 
tracker  to  that  of  the  conventional  tracker  is,  for  bed  <<  n<’  n 


which,  when  jiuJ  is  small,  can  be  further  approximated  as 


Therefore,  since  this  analysis  assumes  that  bed  <<  the  adaptive  filter  bearing 
variance  is  0. 69  times  as  large  (or  1. 6  dB  better)  than  the  split  beam  tracker. 

It  must  be  reiterated  that  this  analysis  Includes  assumptions  whose  validity 
in  a  practical  adaptive  tracker  configuration  have  not  yet  been  verified  by  simula¬ 
tion.  During  the  next  phase  of  this  study,  simulations  will  be  performed,  and  the 
analysis  refined,  if  necessary. 


3. 0  CONCLUSIONS  AND  RELATIONSHIP  TO  PHASE  3 


The  feasibility  of  an  adaptive  filter  aa  a  generic  split  beam  bearing  tracker  has 
been  supported  by  detailed  investigations  Into  certain  of  Its  properties.  In  particular. 
It  haa  been  shown  that  the  static  bearing  estimation  and  dynamic  tracking  perform¬ 
ance  are  comparable  to  or  better  than  conventional  systems,  without  having  to  know 
the  a  priori  spectral  and  temporal  information  needed  to  properly  design  a  fixed 
parameter  device.  The  advantages  of  the  adaptive  filter  tracker  are  most  apparently 
In  the  realm  of  real  world  sensitivity.  In  an  Ideal  environment,  conventional  systems 
are  satisfactory.  Yet  In  actual  applications  they  sometimes  produce  disappointing 
performance.  The  adaptive  tracker  may  provide  the  additional  environmental 
estimates  required  to  perform  closer  to  the  analytical  predictions.  Thus  the  third 
phase  of  the  study  will  emphasize  the  use  of  taped  sea  data  to  drive  simulations 
under  as  many  situations  as  possible.  Analyses  will  stress  the  Inclusion  of 
adverse  environmental  effects  aa  extensions  of  the  analyses  already  performed. 

For  example,  higher  order  dynamics,  multipaths,  and  array  effects  such  as  beam 
steering  and  preformed  beam  handover  will  be  examined. 

A  key  conclusion  of  the  study  of  broadband  plane  wave  Interference  is  the 
bias  effect  produced  on  narrowband  tracker  bearing  estimates.  Broadband  compo¬ 
nents  from  different  directions  within  the  narrowband  of  the  target  energy  do  not 
impact  on  the  variance  of  the  estimate  anywhere  near  a a  significantly  as  they  do 
on  the  mean  value.  This  effect  may  be  the  limitation  In  performance  for  many 
current  narrowband  systems.  The  adaptive  tracker,  which  performs  both  narrow- 
band  and  broadband  tracking  simultaneously  from  the  same  set  of  weight  values, 
may  automatically  resolve  this  problem.  Since  the  data  to  be  employed  Initially 
Is  from  a  towed  line  array  this  effect  will  be  studied  early  In  the  third  phase. 

At  the  conclusion  of  Phase  3  sufficient  real  world  examination  will  have  been 
performed  to  decide  whether  or  not  a  subsequent  sea  test  of  the  resulting  adaptive 
filter  bearing  tracker  would  be  warranted. 
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APPENDIX  A 


Appendix  A 

Bearing  Estimation  Performance  with  Stationary  Narrowband  Signals 


During  the  drat  phase  of  this  study,  the  time  domain  adaptive  tracker  con¬ 
figuration  shown  In  Figure  A-l  was  selected  far  split  beam  bearing  estimation 
with  narrowband  signals.  The  inputs  for  the  narrowband  bearing  estimation 
problem  can  be  modeled  aa 


)->  <t-T) 

d(t)  - 

u 

°s  6 

♦  nd<t> 

H(t) 

x(t)  - 

°s  e 

*  nx(t) 

where  n  ,  and  n  are  independent  aero  mean,  white  gausaian  random  processes 
a  x 

with  power  <r“. 


If  the  target  la  in  the  same  horizontal  plane  as  the  split  array  phase 
centers,  then 


▼  *  £  sin  6  (A- 2) 

with  d  the  distance  between  phase  center,  c  the  speed  of  sound,  and  6 
the  bearing  of  the  target  relative  to  a  normal  line  to  a  vertical  plane 
through  both  phase  centers. 

It  was  shown  during  Phase  1  that  the  steady  state  time  domain  weight 
vector  for  these  inputs  is  a  sampled  sinusoid  with  a  phase  shift.  *  •  -0T. 
To  extract  this  phase  angle,  the  weight  vector  is  Fast  Fourier  Trans¬ 
formed,  and  the  bin  with  the  largest  magnitude  chosen  as  that  containing 


A-l 


Figure  A-l .  Split  Beam  Adaptive  Bearing  Estimator  Structure 


the  signal.  The  phase  of  this  bin,  which  is  corrupted  by  background 
and  algorithm  noise,  is  the  estimate  of  the  phase  angle. 

In  order  to  convert  this  phase  estimate.  4,  to  i  bearing  estimate, 

8.  it  is  necessary  to  know  the  frequency  of  the  signal,  which  is  not 
known  a  priori  and  must  be  estimated.  This  could  be  done  by  inter¬ 
polating  between  the  discrete  FFT  bins  to  determine  the  frequency  of  the 
sinusoid.  However,  it  is  assumed  here  that  the  FFT  bins  are  sufficiently 
closely  spaced  that  the  center  frequency  of  the  bin  with  the  largest 
magnitude  is  taken  as  the  frequency  estimate.  wQ.  The  bearing  estimate 
is  then  given  by 

®  *  iin’l(d^S)  (A‘3) 

The  analysis  that  follows  utilizes  the  frequency  domain  model  for  the 
adaptive  filter  shown  in  Figure  A-2.  During  the  first  phase  of  this  study, 
it  was  shown  that  this  frequency  domain  configuration  behaved  equiva¬ 
lently  to  the  time  domain  filter  for  broadband  inputs,  provided  the  FFT 
time  window  is  large.  The  behavior  was  essentially  the  same  both  in 
transient  and  steady  states. 

Mean  and  Variance  of  the  Frequency  Domain  Weights 

The  primary  advantage  of  the  frequency  domain  model  for  the 
narrowband  tracker  shown  in  Figure  A-2  is  that  It  is  possible  to  compute 
the  mean  and  variance  of  the  weight  vector  with  signal  present.  This  is 
not  possible  in  general  for  the  time  domain  realization  of  the  adaptive 
tracker.  In  Appendix  III  of  reference  [1],  the  mean  and  variance  of  the 
frequency  domain  weights  were  determined  for  a  broadband  input  signal. 
Certain  results  of  that  appendix  can  be  used  here  to  address  the  narrow- 


band  case  under  the  assumptions  that  (a)  the  input  sequences  are 
gaussian,  (b)  the  two  split  array  noise  outputs  are  uncorrelated,  and 
(c)  both  sequences  are  uncorrelated  in  time.  These  assumptions  certainly 
apply  to  the  model  given  in(A-l).  The  reference  also  made  the  assumption 
that  the  inputs  were  zero  mean,  which  is  not  true  here  (the  inputs  have  a 
time  varying,  non-zero  mean).  However,  the  results  used  do  not  require 
this  assumption. 

Using  the  Input  model  given  in  (A-l).  The  reference  also  made  the  assumption 
of  Figure  A-2,  the  Inputs  to  the  single  complex  tap  filter  in  the  kth  FFT 
bin  are 
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or 
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Substituting  gives 


E[W(n+l)1  *  |G(k,w0)r<rs 
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The  steady  state  mean  weight  is  then 
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Again,  using  Appendix  III  of  reference  (1),  with  the  initial  weight 
zero,  the  mean  square  weight  is 

n  n  m- 1 

El  W(n*l)|2)  *  m2  ^  EF(m)i2*2»i2  ^  ^  EtF(m)F*(q)l 

- h*0 -  ms°  c**°  _  iA-10i 

with 
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It  is  shown  in  [1]  that. 

n 
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Using  the  gaussian  assumption,  it  is  easily  shown  that 

Et  d(m)X  (m)|  2]  *  1  G(K.w0>l4  ♦  2M  'G(k.-0)|  2  a  *  *  n*  iA-13) 
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Substituting  this  in  the  first  sum  of  lA-ll)  and  closing  the  sum. 
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Returning  to  reference  (11 
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The  first  fector  Is  evelUble  using  tA-14),  and  It  is  readily  shown  that 
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Finally,  using  iA-14>  from  shove. 
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Substituting  tA-14),  iA-17),  iA-18)  and  v*-19>  I*®  *lvM 
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Closing  the  aumt  gives 


B(F(«)P*(q)J 


To  obtain  the  steady  state  weight,  the  limit  of  (A-15)  and  (A-21)  Is  taken  as 


and  the  results  summed.  L'nder  the  <x>ndltlons,  given  In  (A-24> 


E[|W(.)H 


To  obtain  the  variance  of  the  weight .  subtract  the  magnitude  squared  of 
the  steady  state  mean  weight,  with  the  result  that 

( «•.->  *  •!  »(  4  •*!&*.  .J. m* 

These  statistics  of  the  weight  vector  allow  numerical  evaluation  of  the 
bearing  estimation  bias  and  variance. 

It  is  important  to  note  the  convergence  conditions  for  these  results. 

The  steady  state  mean  and  variance,  (A-9)  and  (A-23),  follow  from  (A- 8)  andA-22). 
respectively,  if 

1  -u(|G(k..0)|2os2  -  Mcn2)  <  1 

and 

l  -  2u(lG(k.-0)i  2a#2  ♦  M?n2)  ♦  u2( |G(k,*0)|  4  c#4  ♦  4|G<k,«-o)|  2Ms(2an2 ♦M2^4)  <1 

Defining  the  signal -to- noise  ratio  in  the  bin  as 

KKk-  JlV* 
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these  requirements  lead  to  the  following  sufficient  conditions; 

0<  uMJn2(tM)«l  iA-24) 
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and 

U  M<^2  (Y2  ♦  4y  ♦  2)  <  2  (Y  ♦  1) 

These  conditions  will  be  important  in  evaluating  the  variance  of  the  bearing 
estimate  in  the  section  to  follow. 


BIAS  AND  VARIANCE  OF  THE  BEARING  ESTIMATE 

The  estimstion  process  consists  of  extrseting  the  phase  from  the  FFT 
bin  with  the  largest  magnitude,  so  it  is  necessary  to  obtain  the  density  of 
the  phase  in  each  FFT  bin.  Because  the  input  statistics  are  gaussian  and 
a  large  number  of  iterations  of  the  weight  recursion  equation  can  be 
assumed .  the  frequency  domain  weights  may  be  regarded  as  complex 
gaussian  random  variables.  Therefore,  the  density  of  the  phase  of  the  k**1 
frequency  domain  weight  can  be  obtained  by  transforming  the  two  dimen¬ 
sional  gaussian  density  of  the  real  and  imaginary  parts  of  Wk  to  polar 

f  21 

coordinates.  This  gives 
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with 
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The  phase  density  is  obtained  by  integrating  with  respect 


to 
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The  variance  of  the  phaae  can  be  obtained  by  numerical  integration  of  (A-26). 
The  result  is  parametric  in 


which  for  the  adaptive  tracker  ia  given  by 
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For  a  real  sinusoidal  input,  the  FFT  power  gain  la  no  laager  given  by  (A-27),  and 
the  expreeaiona  far  the  variance  are  considerably  more  complicated.  For 
example,  while  for  the  complex  case, 

M-l  M-l  Ja  (nTs-r)  Ju^a  . TT  (n-m)k  -  -Ja  r 
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the  real  case  gives 
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Fortunately,  when  u»o  is  near  bin  renter,  the  results  given  above  hold  if 
j  G(k,  WQ)j  2  is  replaced  by 

|Gr(k.-0)l  2  “  ll  G<k-0)|  2  (A’32) 

and  the  signal  power  is  appropriately  adjusted. 

Of  interest  are  the  statistics  of  4  rather  than  that  of  the  phase. 


The  bearing  estimate  Is  given  by 


where  1>  Is  the  center  frequency  of  the  bin  with  the  largest  magnitude, 
o 

Assuming  that  sufficient  integration  la  used  to  assure  selection  of  the  bin 


cloeest  to  the  signal,  call  it  bin  k. 


e  *  sin 
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with  * k  the  phase  of  the  k  bin.  Then 
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Var  l?l  *  E  l$2l  -  E2  [•). 


While  these  Integrals  cannot  be  evaluated  explicitly,  they  are  readily 
evaluated  numerically  on  the  computer. 

Some  care  is  required  in  interpreting  the  bearing  statistics  determined 
by  this  method  due  to  the  modulo  2»  assumption  inherent  in  the  derivation 
of  the  phase  density.  The  bearing  estimate  given  by 


8  *  sin 


*  /CMJl  t  \ 

\2»kd  #xy 


has  considered  *k  to  be  on  t-».»  ).  that  is.  modulo  2  *.  so  assuming  the 
principal  value  of  the  arc  sine .  the  bearing  estimate  is  restricted  to 
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The  statistics  of  th«  boaring  aatimate  aa  derived  above  do  not ,  therefore . 
reflect  the  true  bearing  error  statistics  when  the  phase  makes  frequent 
excursions  of  more  than  *  radians.  This  is  similar  to  the  consideration  of 
cycle  skipping  in  the  analysis  of  phase  errors  in  phase  locked  loops 

The  phase  density.  (A- 26),  is  shown  In  Figure  A-3  aa  a  function  of  the 

2  2 
parameter  s  .  It  can  be  seen  that  for  the  parameter  8  >1.  the  density  is 

essentially  zero  near  *  >it,  so  such  excursions  are  unlikely.  Under  this 

condition,  the  bearing  estimate  statistics  evaluated  by  the  method  given 

above  should  be  valid.  Figure  A-4  shows  the  standard  deviation  of  *as  a  func- 

O  *> 

tion  3~ ,  computed  numerically  from  (A-26).  For  B~  •  1,  the  phase  standard 
deviation  la  48°,  If  9  is  small,  this  give  a  bearing  error  standard  deviation  of 
approximately  (48  c/u>od).  This  is  a  nominal  upper  bound  on  the  range  over 
which  the  results  of  this  section  are  valid. 

Figure  A- 5  shows  the  predicted  rms  bearing  error  for  a  particular  real 
adaptive  tracker,  computed  as  described  above.  The  split  array  has  150  ft 
between  phase  centers,  and  the  signal  la  a  100  Hz  real  sinusoid  embedded  in  white 
gausslan  noise.  The  FFT  has  128  points  and  a  resolution  of  2  Hz.  so  the 
signal  is  centered  in  the  50th  bin.  The  figure  shows  the  rms  bearing  error 
as  a  function  of  the  signal- to- noise  ratio  at  the  FFT  input  for  various  values 
of  the  weight  feedback  coefficient .  u  .  In  keeping  with  the  comments  of 
the  preceding  paragraph,  the  curve  is  valid  for  rms  bearing  error  of 
less  than  2.54°.  Further,  the  curves  are  only  plotted  over  .the  range  of 

t 

signal- to- noise  ratio  for  which  the  stability  conditions  of  the  preceding 
section  are  met  accounting  for  the  endpoints  on  the  curves  for  *  *  2  *  * . 
2"12,  and  2~13.  Computer  simulations  over  10  ensembles  are  also  shown  on 
the  figure,  supporting  the  analytical  predictions. 
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Figure  A-J.  Phase  Density  for  Various  VaJues  of 


Figure  A4  Standard  Deviation  of  Phase  for  Density  of  F^rr*  3 
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Figure  A-5.  Predicted  RMS  Bearing  Error 
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AN  APPROXIMATION  TO  THE  BEARING  ESTIMATE  STATISTICS 

In  the  special  case  of  high  signal-to-noiae  ratio  and  small  bat  ring  angle, 
U  is  possible  to  develop  a  simple  approximation  to  the  mean  and  variance  of 
the  bearing  estimate.  At  small  bearing  angles 


a  ■  sin 


'cMTg  "]  cMTb 
2  *  kd  *k  *  2t  kd 


lA-54) 


so  the  bearing  statistics  can  be  readily  related  to  those  of  the  phase  angle 


The  phase,  #k>  is  given  by 


♦k  Stan 


lA-35) 


where  is  the  complex  weight  in  the  Kin  FFT  bin.  From  above.  Afc  and 
are  uncorrelated  gauasian  random  variables  with  mean 


1°  <*~0),2°e2  d 

- 2_ — -» - -  gin  t  —  «in  9 ) 

a  it  •  2  ,ii  *  OC 


E  [A  )  *  - — f - 

^  !0(k.-0)|2a,2^n 
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|G  (k.-  )12  3  2 


E  (B.  ]  »  - “-5 - f - j  C°l  (*o78ln  9> 

k  |G(k.0)!2c#2^n2 


Var(AkJ  *  Var  [B^j  -  (1/2)  Var  [WjJ 
where  VnrtW^]  is  given  in  (A-23). 
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It  Is  possible  to  epproximste  the  expectation  of  a  function  of  two  random 
variables,  g(x.y),  given  the  Joint  density  function  of  x  and  y,  Pxy(**y) 
by  expanding  g(x.y)  In  a  Taylor  series. 


g(x.y)  ■  g(a.b)  ♦  [|f  <*  _  •>  *  j*<y  ‘  b>j 

♦  2  (x  -  a)  (y  -  b)  ♦  -***■  (y  -  s)2l  < 

*  y  jy  J 


(x  -  a)4 


where 


(x.y> 


htn"r  «xr  3xn'ra  yr 
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If  a  »»7  and  b  ■»  .  the  means  of  x  and  y.  the  Joint  density  p  (x.y)  is  con- 
x  y 

centrated  near  end  g(x.y)  is  smooth  in  the  vicinity  of  P7x.»]y)*  then 

E(g(x.y)l  can  be  approximated  by  retaining  second  order  terms  in  lA-39) 


giving* 

Elg(x.y)!  .  »(ns.ny)  ♦  *  ♦  2uJfy  ♦  cy2  l -fj 


VV-40) 


with  ■.  *  E  ICx  -  n  )(y  -n  )].  SimUarly .  by  expanding  g  (x  .ylln  a 

*y  *  y  . 

Taylor  series  and  retaining  only  second  order  terms,  it  can  be  shown  that 


2  at  H  || 


°y  * 


ax  ay  xy 


(A-41) 
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Th#  cm*  of  interest  her*  Is  g(x.y)  -  tsn'1  with  y  ■  Afc.  x  •  B^,  so 

»>y  ■  Ely)  ■  CjSin  (-2-sin  s).  nx  ■  E(x)  •  CjCoa^^  sin  fl)  (A-42) 


2  2 
G(k,  wQ)  o& 


2  2  'J'*»  “o 

o  .  ff  -  VsrtAJ.  u  -  0.  and  C  -  - 1 - - - 5 - ;>  . 

y  ^  1  |c«-  “of*  *.  ♦  «% 

To  evaluate  these  approximations,  w*  need  th*  derivatives  ig/}x.  jg/  ay. 
*  2g/  >x*.  d“g/dy”  are  needed. 

Car*  must  be  taken  in  th*  interpretation  of  g(x.y)  if  the  result  is  to 
be  valid  for  Sc  l  -  *.  »}.  This  is  due  to  the  fact  that 


a  «  tan  l(a) 


is  generally  Msumed  to  return  ac  [>  j  .  >] .  Th*  estimator,  on  the  other 
hand,  returns  •«(-*.*)  by  recognizing  th*  signs  of  x  and  y.  Therefore, 
g(x.y)  is  properly  defined  as 


g(x.y)  •  nan*1  f£ 


where  *tan(*)  Includes  th*  quadrant  resolution  using  th*  signs  of  x  and  y. 
Using  ny ,  and  cy  from  vA -42),  the  derivatives  in  tA-40)  and  (A-41)  are 


>1  .  . 


C  jStnx 


(CJ  sin  x  ♦  Cj  cos*  x ) 


H  . 

>y 


C  jCosx 


(Cj  sin  x  ♦  Cj  cos'x) 


*  g  2(C  jCOs  x)C  jSin  x 

2  "  ”1 - 1 - 2 - 2 — 

»x*  (C*  sin*X  ♦  Cj  sin'  X) 

,2 

g  sin  2  \  d  _ 

ay2  cT  oc 
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sin  2  x 
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Therefore, 


«'V  -  -  ‘  [|g£]  *  i  £  ‘  -  v*  — ]-'4 


sin«  lfi-44) 


<*«*  2  av'  2  WX~ 

Var  (#kJ  ■  — y  •*»  *  ♦  -y  cos  A  *  -y 

C,  C ,  C. 


(A -45) 


Using  the  small  baarlng  angle  approximation.  tA-34),  we  get 

El#l'  iTB*-ElM 


and.  using  the  definition  of  9  from  |A-29), 


«5T  1 


v„  s;  -L 

(2»k)  d*  29 


(A -46) 


(2>3  ♦  4v2  ♦  4>  ♦  1) 


-  c  _ (gy  ♦  4y  »  n  ♦  n _ 

(2«k)5  d2  y2[2(y  ♦  1)  -«M  on2(Y^  «*♦  *>]<?  ♦  1>TJ 

This  approximation  is  compared  with  the  numerically  evaluated  r.m.s. 
bearing  error  for  a  particular  real  adaptive  filter  in  Figures  A-6  through  A-8. 

The  tracker  uses  a  128  point  FFT  with  two  Hx  resolution.  The  100  H*  real  stnu- 

|A. 

aoidal  signal  embedded  in  white  noise  is  centered  in  the  50  FFT  bin.  The 
results  are  given  for  three  values  of  p,  and,  as  before,  the  curves  are  only  plot¬ 
ted  over  the  range  of  SNR  for  which  the  filter  Is  stable.  The  target  bearing  Is 
aero  degrees.  As  can  be  seen,  the  approximation,  (A-46),  gives  excellent 
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agreement  with  the  numerical  computation  over  the  range  for  which  the  analyala 
la  valid,  that  la,  for  atandard  deviation  leaa  than  2. 54  degreea.  Thla  auggeata 
that  for  a  mall  bearing  anglea,  the  numerical  computation  can  be  avoided. 
SIMULATIONS  OF  THE  NARROWBAND  ADAPTIVE  TRACKER 

Slmulatlooa  of  the  narrowband  caae  were  performed  for  two  array 
configurations,  each  with  a  different  signal  frequency  embedded  in  white  Gaussian 
noise.  A  vernier  FFT  with  a  resolution  of  2  Hz  is  positioned  around  the  location 
of  the  narrowband  signal,  with  the  signal  centered  In  one  of  the  frequency  bins. 
The  phase  of  that  bin  la  extracted  and  converted  to  a  bearing  estimate.  The  com¬ 
puter  calculates  the  mean  and  variance  of  the  estimate  using  a  sliding  window, 
and  generates  a  plot  of  the  estimate  vs.  time. 

The  simulation  model  la  a  PDP-11/70  computer  implementation  of  the 
tracker,  driven  by  pseudo-random  noise  and  signal  generators  designed  to  simu¬ 
late  environments  of  interest,  as  shown  In  Figure  A-9.  Three  Independent 
pseudo-random  gaussian  noiae  (PRGN1  generators  provide  independent  background 
noise  for  the  simulated  half  array  outputs  and  the  broadband  signal.  These  are 
each  passed  through  a  finite  impulse  response  filter  with  selectable  response, 
allowing  simulations  of  various  signal  and  noiae  spectra.  The  Inter-array  delay 
can  be  injected  either  as  a  fixed  value,  or  can  vary  with  time.  In  addition,  a 
sinusoidal  signal  component  can  be  added  to  the  simulated  array  output,  with  the 
appropriate  phase  shift  between  array  halves.  The  model  also  includes  the 
capability  to  plot  the  resulting  estimates,  and  to  compute  second  order  statistics 
using  s  sliding  window  technique. 

The  first  array  simulated  has  150  feet  between  phase  centers,  and  a  100  Hz 
line  emanating  from  a  target  at  0°  bearing  relative  to  broadside.  The  vernier 
16  point  FFT  ta  positioned  such  that  the  line  is  centered  in  FFT  bin  number  four. 
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Figure  A-9.  Simulation  Model  For  Adaptive  Tracker  Input* 


Figures  A- 10  through  A- 12  show  the  bearing  estimate  vs  time  for  a  feedback 
-12 

coefficient,  p,  of  2  and  Input  slgnal-to-nolse  ratios  of  9,  -1,  and  -11  dB. 

The  Input  SNR  Is  measured  In  the  entire  FFT  Input  band,  and  the  selection  of  p 
yields  a  time  constant  of  128  seconds  for  the  adaptive  filter.  The  simulations  are 
approximately  8  time  constants  In  duration.  Figures  A-13  through  A- 16  show  the 
same  data  for  p  =  2_1*  and  SNRs  of  9,  -1,  -11,  and  -21  dB.  In  this  case,  the 
time  constant  la  512  seconds,  so  the  runs  are  about  2  time  constants  long.  Fig¬ 
ure  A- 17  compares  the  simulation  results  with  theoretically  predicted  values.  In 
all  cases,  the  statistics  were  computed  over  the  270  most  recent  estimates. 

The  other  array  simulated  has  75  feet  between  phase  centers  and  a  600  Hz 
signal  from  a  target  at  broadside.  As  above,  the  16  point  FFT  Is  positioned  such 
that  the  signal  Is  bin-centered  In  frequency  bin  number  four.  Figures  A-18 
through  A- 20  show  the  bearing  estimate  vs  time  for  p  *  2~10  and  SNRs  of  10,  0, 
and  -10  dB.  The  time  constant  Is  32  seconds,  and  the  simulations  are  about  32 

time  constants  In  duration.  Figures  A-21  through  A-23  give  the  results  for 
-12 

M  *  2  ,  and  for  the  same  SNRs,  corresponding  to  a  128  second  time  constant 

and  a  run  duration  of  8  time  constants.  The  results  for  g  »  2"1*,  with  SNR's  of 
10,  0,  -10,  and  -20  dB,  respectively,  are  shown  In  Figures  A-24  through  A- 27. 
These  have  a  time  constant  of  512  seconds  and  simulation  length  of  approximately 
4  time  constants.  The  simulation  results  are  compared  to  the  theoretical  predic¬ 
tions  (computed  as  described  earlier)  in  Figure  A-28,  where  statistics  were 
computed  over  the  200  most  recent  estimates. 

Finally,  simulation  runs  were  made  using  the  tracker  configuration  used  in 
the  predictions  of  Figure  A-5.  Since  the  analysis  has  a?*um?d  that  the  FFT  bln 
containing  the  signal  Is  selected  for  phase  extraction,  it  was  felt  that  the  best 
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Fifuit  A- 13.  Bearing  Eatimatc  Vcr*u»  Time 
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Figure  A  19.  Bearing  Estimate  n  Time 
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Figure  A-21.  Bearing  Emuiwic  Versus  Time 


Figure  A- 24.  Bearing  Estimate  Vtmu  Time 
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check  on  the  analysis  would  be  to  force  the  simulation  to  select  that  bin.  There¬ 
fore,  for  the  simulation  points  shown  on  Figure  A-5,  the  tracker  was  forced  to 
always  choose  FFT  bln  50.  containing  the  100  Ha  signal.  This  may  not  be  as 
unreasonable  from  an  operational  point  of  view  aa  it  first  sounds.  If  a 
conventional  FFT  detector  with  gain  exceeding  that  of  the  tracker  FFT  were  used 
to  determine  line  frequency,  the  tracker  could  operate  in  the  correct  bln  at  low 
SNR.  The  simulation  results  shown  in  Figure  A-5  are  ensemble  averaged  over 
10  tracker  runs,  aa  well  as  time  averaged  over  190  estimates  within  each  run. 


A -48 


Appendix  B 


Comparison  of  NB  Adaptive  Tracker  Performance  with  the 
Cramsr-ftao  Lower  Bound 


The  Cramer-Reo  Lower  Bound  (CRLB)  is  a  lower  bound  on  the  variance 
of  any  eatimate  of  a  parameter  given  a  particular  aet  of  obaervatlona.  For 
the  eatlmation  of  multiple  parametera .  or  the  eatlmation  of  a  particular 
parameter  when  othera  are  unknown  (and.  hence,  muat  be  eatlmated)  the 
CRLB  can  be  determined  using  Fisher's  information  matrix. ^  Let  Z  be  a 
vector  of  observations  and  A  be  a  vector  of  parameters  to  be  estimated . 
where 


A  *  [Aj.Aj.  Ag] 


(B-l) 


Fisher's  information  matrix  Is  a  K  X  K  matrix  given  by 


J  ■  -  El  <2at lln  Px/a  (Z/A)))l 
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where 


Px>#(Z/A)  ■  conditional  density  of  Z  given  A 


r 2i  _iiiT 

1“.  *  »\J 
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Any  unbUMd  estimate  of  A{  satisfies 


Var  [Ajd)  -  AtJ  *  (B-4) 

where  is  th«  li**1  element  of  J  * . 

Lot  the  2N  dimensional  observation  vector  for  ths  narrowband  bearing 
estimation  problem  in  the  kth  frequency  bin  be  given  by 

^  «  [Dk(l).Xk(l).  - —  Dr(N).  Xk(n)]T  (»-») 

where  Dk  and  Xk  are  the  FFT  outputs  In  the  bin  closet  to  the  signal  as 
defined  in  <A-5>.  Assume  that  the  signal  amplitude,  o 8.  frequency,  wo,  and 
bearing,  9,  are  unknown  and  to  be  estimated.  Then 

A  *  tog.  <*>„.  #1T  l®-6* 

Using  the  mean  and  variance  of  Dk  and  Xk  from  (A-7),  and  the  gausslan  input 
assumption .  the  conditional  density  of  Z  given  A  is 
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where  ♦  denotes  complex  conjugate  transpose  sod 

[  -J"0  |  •  -J"o  c  * 

k(A)  ■  E[Z  A)  *  o§S(k.^)  Le  °  c  .1, - « 


B-2 


Tharafora , 
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Now.  obMrvt  that 
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Not*  that  this  matrix  ia  Harmltian  and  danota  its  l)in  alamant  as  3^. 
Tha  alamants  of  J  can  ba  found  by  considarlnp  tha  darlvativas 
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Substitution  of  th«s«  dsrivativs*  to  yisld  lu*  sJsinsnts  a  ^  jivss 
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The  term  of  Interest  Is  J  ,  the  i»3,  )«3  element  of  J  ,  since 


Vsr  (Aj  (Z)  -  AjJ  «  Vsr  [6  <Z)  ■  «]  1 


This  element  is  given  by 


where .  substituting  from  (B-17),  after  some  algebra 
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Therefore 


where 


M  3. 
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Then  J  ,  which  is  the  CRLfi  for  the  boaring  •stlmata,  la 


.33 
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This  can  b«  expressed  In  terms  of  tha  aignai  to  noiaa  ratio  in  tha  bin 


S2a  2 

?  a.L 


|  S(k.MQ)|  2ea2 


M  o. 
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so  that 


Var  [9  -  9]e 
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yN  *q2  -^2  co§29 
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1  *  c  ,ine  1  *  (M-'i)f  c  §in  e]]  (B“26) 

a  •*  a 


It  is  interesting  to  compare  this  result  with  tha  high  SNR  approximation 
tor  tha  variance  of  tha  adaptive  tracker  bearing  estimate.  iA-46)  developed  In 
the  preceding  section.  For  small  bearing  angles  and  high  SNR,  the  variance 

la 


M2T  2 

Var  t «]  • - *— 

(2  «k) 


(B-27) 


where,  under  the  same  assumptions.  (B-26)  reduces  to 


Var  [9  -  91  *  — 


(B-28) 


Sine*  MTg/k  is  the  estimate  of  uQ,  it  can  b«  seen  that  thaie  is  a  correspondence 
between  the  time  constant  of  the  adaptive  filter.  In  its.  .tons,  and  the  total 
number  of  observations  of  the  FFT  output  in  the  CRLB  as  would  be  expected. 
Just  as  in  the  case  of  the  comparison  of  the  adaptive  tracker  with  the 
CRLB  for  broadband  Inputs. ^  it  is  necessary  to  resolve  the  relationship 
between  N  and  t..  The  method  applied  in  the  broadband  analysis,  specifi¬ 
cally.  the  incorporation  of  an  exponentially  decay  term  In  the  large  SNR 
approximation,  also  works  well  here. 

The  mean  weight,  given  a  *ero  initial  weight,  can  be  written  as  a 

function  of  time  as 

ElWk(n)l  -  EtWk)[l  -  e  **] 

where  n  is  the  number  of  filter  iterations  and  ElWjJ  is  the  steady  state 
mean  weight .  The  variance  has  approximately  half  the  time  constant .  so 
we  can  assume  that  after  one  or  so  time  constants .  the  transient  response 
of  the  variance  may  be  neglected.  Then  Incorporating  this  change  in  the 
derivation  of  (A-*®)  gives 

mH  2  2  .  r  -*»/t  ii 2 

Vw  ItcnH - 4  2,  ‘  ll  -  .  ']  (B-JO, 

(2*k)Z  d2  ^iL  J 

The  parameter  n  corresponds  exactly  to  the  number  of  FFT  outputs 
processed  by  the  filter,  that  is,  to  N  in  the  CRLB.  Then 


It  is  assumed  here  that  the  FFT  resolution  is  sufficiently  fine  that 

m  m  2*  k/MT  .  This  ratio  is  plotted  as  a  function  of  the  ratio  N/  r.  in 
os  .1 

Figure  B-l.  It  con  be  seen  that  the  minimum  occurs  at  N/  r(  «  1,25,  and  that 
the  ratio  is  1.228  at  the  minimum.  Therefore,  at  high  SNR.  the  adaptive 
tracker  is  degraded  approximately  0.45  dB  from  the  CRLB.  The  relative 
performance  at  low  SNR  must  be  determined  numerically  as  described  earlier. 
Whether  or  not  the  ratio  SI  Tj  «  1.25  yields  the  best  relative  performance 
at  low  signal  to  noise  ratios  and  at  other  bearing  angles  must  be  determined 
numerically . 

Figures  B-2  through  B-4  compare  the  performance  of  a  particular  real 
input  narrowband  adaptive  tracker  with  the  CRLB  for  a  100  Hz  real  sinusoidal 
signal  from  a  target  at  0°.  The  FFT  size  is  128  points  and  the  resolution  is 
2  Hz.  The  distance  between  the  split  array  phase  centers  is  150  feet.  Three 
values  of  the  feedback  coefficient,  n,  are  shown,  and  the  curves  only  plotted  in 
the  stable  range  of  the  filter.  The  dashed  portions  of  the  predicted  rms  error 
curves  represent  the  range  in  which  the  numerical  analysis  is  not  valid. 

It  can  be  seen  that  above  an  input  SNR  of  -20  dB .  corresponding 
nominally  to  a  0  dB  SNR  in  the  bin.  the  ratio  of  the  predicted  rms  error  to 
the  square  root  of  CRLB  is  on  the  order  of  1 . 15  to  1 . 5  if  u  is  small  enough. 
This  agrees  well  with  the  ratio  of  1.228  developed  by  approximation  above. 
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Fifurf  B-3.  Companion  of  Predicted 


Error  with  Cramer-RAO  Lower  Bound 


•  IMMI 


Appendix  C 


Mean  and  Variance  of  the  Adaptive  Tracker  Weights  for 
Simultaneous  Broadband  and  Narrowband  tHrcctional  Source 


The  adaptive  Alter  configuration  of  Figure  C-l  haa  beon  proposed  aa  an 
alternative  to  conventional  split  beam  trackers  as  part  of  the  Adaptive  Tracking 
System  Study  In  this  *PPr°ach.  the  split  array  outpuU  are  provided  as  the 

primary  and  reference  inputs  to  the  adaptive  Alter,  fa  the  case  of  broadband  sig¬ 
nals,  the  bearing  of  the  target  is  determined  by  extracting  the  time  delay  of  the 
target  wavefront  between  array  halves  from  the  time  domain  weight  vector. 

This  is  done  by  locating  the  peak  of  the  weight  vector,  its  location  corresponding 
to  the  delay.  For  narrowband  signals,  the  phase  shift  between  array  halves  Is 

extracted  from  the  Fast  Fourier  Transform  (FFT)  of  the  weight  vector  and  coo¬ 

rd 

verted  to  bearing. 

Previous  results  have  considered  the  mean  and  variance  of  these  estimates 
for  broa<t>and  signal  only^  and  narrowband  signal  only  in  Appendix  A.  It  la  now 
desired  to  consider  (1)  single  targets  with  both  broadband  and  narrowband  compo¬ 
nents,  and  (2)  simultaneous  broadband  and  narrowband  sources  from  different 
directions,  one  the  target  of  interest  and  one  an  Interference.  Analysis  of  the 
adaptive  tracker  performance  In  this  case  requires  determination  of  the  mean 
and  variance  of  the  weight  vector  with  these  Inputs. 

These  statistics  are  developed  in  this  Appendix  by  using  the  frequency 
domain  realisation  of  the  algorithm,  as  done  In  Appendix  A  for  narrowband 


C- 


stgnaia  only. 


Let  the  left  end  right  half  arrey  inputs  to  the  adaptive  tracker  be 


d(t)  ■»  0^^°^  ^  ♦  b(t  -  Tj)  ♦  n^lt) 

Iw* 

x(t)  *  <r  e  ♦  b(t)  ♦  n  (t) 

8  * 


(C-l) 


where  n  .  n. .  and  b  are  Independent,  zero  mean,  white  gauasian  random 
x  d 

processes.  The  delay,  r  .  is  the  delay  of  the  narrowband  signal  between 
array  halves,  given  by 


and  t,  the  delay  associated  with  the  broadband  signal. 


r 

2 


(C-2) 


(C-S) 


where 

d  =  distance  between  array  centers 

c  *  speed  of  sound 

g  »  bearing  angle  of  narrowband  source 
9  •  bearing  angle  of  broadband  source 

It  is  assumed  that  the  targets  are  in  the  same  horlsontal  plane  as  the 
array  centers,  and  that  r2  is  an  integer  multiple  of  the  sample  interval  of 
the  algorithm.  T#.  that  is 


Reference  ( 1)  considered  the  mean  and  variance  of  the  weight  vector 
of  the  adaptive  tracker  by  utilising  a  frequency  domain  model  for  the 
adaptive  filter.  This  has  been  shown  to  behave  equivalently  to  the  time 
domain  adaptive  filter  when  the  FFT  time  window  is  large.  Appendix  A 
extended  these  results  to  the  case  of  a  narrowband  input  signal.  The 
assumptions  necessary  to  use  these  results  here  are:  (a)  gaussian  input 
sequences,  (b)  noise  uncorrelated  between  split  array  outputs,  and  (c) 
frequency  domain  input  sequences  un correlated  in  time.  It  is  shown  in 
below  that  (a)  and  (b)  are  true  and  that  (c)  is  approximately  true  for  large  FFT 
sizes. 

The  inputs  to  the  single  complex  tap  adaptive  filter  in  the  k**1  FFT 
bin  are 


Xk(n) 


m-1 

ito 


nhl)  T#)  e 


and 


(C-5) 


m-1 

Dk(n)  *  £  dtff  -  nil)  Tf] 

1*0 

The  FFT  is  assumed  to  be  performed  with  no  redundancy  and  no  gap 
between  successive  time  records.  These  can  be  written  aa 


Xk(n)  *  G(k,  wQ)  <7a  ♦  Nxk(n)  ♦  Bxk  (n) 


(C-«) 


C-3 
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and 


Dk(n> 


G(k,  u»Q)  <t8* 


-J  <Vi 


♦  Ndk(n) 


♦  B  (n) 


where  Nxk.  Bxk.  N^,  end  B^  ere  the  tranaforma  of  b(t),  nd«  and 
b(t  -  t  j)  t  reapectlvely .  and 


G(k.  uQ) 


•4f<¥k  -^T,)]  ¥  (  Mk  “»T<) 

r  •  ft  _  •% 


*ln[T(2Sk-"oV] 


j  W  nMT 
‘  o  a 
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It  la  aasumed  that  k  **0,  M/2  since  the  reaulta  differ  In  these  bins.  It  was  shown 
In  reference  [1]  that  for  k  "0,  M/2, 

EJN^  -  e;  ^  -  0 

(C-8) 

EINxkNdkJ  “  0 
and 


E(  N*k(ro>Nxk(p>I 
El  Sdk(m,Ndk(p’l 


j  o,  p  0  m 
{  Mum2,  p  -  m 

|  o,  p  -  m 
j  p  -  m 


<C-9f 


It  la  alao  easily  shown  that 

EJB^fHM)  -  EIB^fmi)  -  0 

■  bIba("“b*«»1  -  I  Mc^.p-n. 

However, 

Mm1  M-l  2r 

EfB^nnB’fpil  -  I  I  E{b(  ( J  bJfq-pnnT  -t«J 
f-0  q-0 

«  <T.2  I  I  6«£-q>T.+  (m-piMT-r^e'1**  {i~*)k  (C-10J 
£-0  q-0 
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Given  the  range  of  /  and  q,  if  in-pl  >  1,  then  the  expectation  is  zero.  If  m  •  p,  then 


•  2  .|2l  AW 

EfB^lmiB^fp*)  -  <rb  <M-A>e  * 


(C-ll) 


and  if  Im-pl  -  1 

ElB^nnB^p*!  -  (Tb2Ae'JS  (M^)k 


.  2* 


Jh, 


(C-12) 

The  input  samples  to  the  single  tap  filter  in  the  klu  bln  are  not,  therefore,  uncor¬ 
related  aa  assumed  In  Cl3.  However,  the  correlation  coefficient  between  adjacent 
samples  (im-pi  ■  1)  is 


m-\A 


(C-13) 


Therefore,  if  M  >  >  A  ,  the  input  sequence  is  approximately  uncorrelated,  and  we 
the  resulta  of  Appendix  111.  reference  [2]  apply. 

From  (2],  the  mean  weight  in  the  kth  bin  can  be  written  as 

n  »  n  2. 

F.[W(n+l)]  -Ml  E(  D.  (m)X.  (n»)  n  E(l-MX(kjn  (C-14) 

m-0  k»m>l 


Using  the  above  model, 
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and 
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Substituting  in  (C-15)  gives 
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n  ♦! 


lG(k,u/0>  |  <tb  Mf  ♦  ob  > 


.  £  JO(k, u)q) 1 2<rs2 e'^oH  +  (M-A )  obVJ>l  Ak] 
Under  the  condition  that  j  l-*i  |  |G<k, 4oq>|  2Gg2  ♦  M(on2^  <rb2»)J  <  1  , 


(C-16) 


thla  gives  a  steady  atate  mean  weight  of 
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Reference  (1]  also  frivea 


B|F(m)|2  =  E  [|Dk(«n)  XR(n)|  2]  TT  ■([  1  -  M|Xk(p)|  2]  I 

p=tn+l  l  / 

With  considerable  algebra,  it  can  be  shown  that 
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with 


M2  *  |  G(k.«0)|2 
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m=0 

and  if  ;A4;<  1.  thia  goes  to 

»*  ^  E  [|  F  ( n )  |  2]  *  *2  X1 

M=0  4 

To  evaluate  the  second  term  of  (C-18).  the  evaluation  of 
E  lFk<m)  la  required  *l*. 

E  lFk(m)  F*(q)J  *  E  [D*(q)  XR(q)J 
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Therefore 
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Substituting  this  in  ( C-25 )  and  closing  the  sums  gives 
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If  hP^.  <  I  and  A4  j  <  1  .  as  n  —  •  ,  this  goes  to 
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The  mean  square  weight  in  steady  state  ia  therefore 
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To  obtain  the  variance,  subtract  E(W(-)jj  2.  and  obtain 
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2  2 

When  either  *  0  or  12  ”  =  0.  this  checks  with  previous  results. 

In  most  esses  of  interest,  either  the  tracker  will  be  implemented  in 
the  time  domain,  so  that  no  windowing  occurs,  or  M  »A,  so  windowing 
is  negligible.  In  that  case 


M  -  A  «  M 


and 


S(M.  A)  •  0 
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Therefore 
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and  Aj  is  unchanged .  This  can  be  rewritten  a. 
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Note  that  the  ooalne  term  vaniahee  whenever  =  2Ak/M(  which  happen, 
when  the  broadband  and  narrowband  energy  come  from  the  aame  direction 
and  the  narrowband  signal  ia  bin  centered. 


Summary 

Expresalons  (or  the  mean  and  variance  o l  the  complex  weight  of  the  frequency 
domain  adaptive  tracker  In  the  presence  of  simultaneous  broadband  and  narrow- 
band  directional  sources  at  arbitrary  bearing  angles  have  been  developed. 

These  results  allow  determination  of: 

( 1)  Performance  of  narrowband  adaptive  tracker  algorithm  for 
targets  with  broadband  and  narrowband  components  by  the 
method  of  Appendix  A. 

( 2)  Effect  of  broadband  interference  on  narrowband  tracker  by  the 
same  method. 

( 3)  Performance  of  the  broadband  tracker  algorithm  for  targets 
with  narrowband  components  using  the  method  of  (11. 

(4)  Effects  of  narrowband  interference  on  the  broadband  tracker  by 
the  same  methods. 
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Pifurc  C  l.  Adaptive  Tracker  Structure 


During  Phase  1  of  the  Adaptive  Tracking  System  Study,  the  adaptive  filter 
configuration  of  Figure  C-l  has  been  proposed  as  a  means  of  determining  the 
bearing  of  a  narrowband  source.  The  outputs  of  the  two  split  arrays  provide  the 
reference  and  primary  Inputs  of  the  adaptive  filter.  When  only  the  narrowband 
target  and  the  uncorrelated  hydrophone  noise  are  present,  the  mean  weight  vector 
converges  to  a  sampled  sinusoid  at  the  same  frequency  as  the  narrowband  source 
with  a  phase  shift 

«  -  ^  |  sin  *  (D-l) 


where  -0  is  the  source  frequency.  9a  is  the  source  bearing  angle,  d  the  distance 
between  array  phase  centers,  and  c  the  speed  of  sound.  It  is  assumed  that  the 
target  is  in  the  same  horizontal  plane  as  the  array  phase  centers.  To  extract 
this  phase  angle,  the  weight  vector  is  Fast  Fourier  Transformed,  and  the  bin  with 
the  largest  magnitude  chosen  as  that  containing  the  signal.  The  phase  of  this  bin, 
which  is  corrupted  by  background  and  algorithm  noise,  is  the  estimate  of  the 
phase  angle. 

.  A 

In  order  to  convert  this  phase  estimate,  o  ,  to  a  bearing  estimate,  #,  it  is 
necessary  to  know  the  frequency  of  the  signal,  which  is  not  known  a  priori  and 
must  be  estimated.  This  could  be  done  by  interpolating  between  the  discrete  FFT 
bins  to  determine  the  frequency  of  the  sinusoid.  However,  it  is  assumed  here 
that  the  FFT  bins  arc  sufficiently  closely  spaced  that  the  center  frequency  of  the 
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bln  with  the  largest  magnitude  la  taken  as  the  frequency  estimate,  -»Q.  The 
bearing  estimate  Is  then  given  by 


(D-*> 


Appendix  A  developed  the  mean  and  variance  of  the  bearing  estimate  of  such 
a  tracker  when  the  narrowband  source  Is  the  only  directional  component  In  the 
noise  field.  Considered  herein  is  the  case  when  a  broadband  directional  source  is 
also  present.  When  located  at  an  angle  different  bom  the  narrowband  source 
this  may  represent  a  broadband  Interference,  while  the  broadband  source  at  the 
same  angle  as  the  narrowband  source  models  a  single  target  with  components  of 
both  types.  In  Appendix  A,  the  mean  and  variance  of  the  bearing  estimate  for 
the  narrowband  adaptive  tracker  with  a  sinusoidal  signal  was  determined  as 
follows: 


(1)  Derive  the  mean  and  standard  deviation  of  the  FFT  of  the  time 
domain  weight  vector. 

(2)  Based  upon  the  large  number  of  iterations  required  in  the  adaptive 
algorithm,  assume  that  the  FFT  of  the  weights  are  complex  gaussian 
random  variables. 


(3)  Transform  the  weights  to  polar  coordinates,  then  integrate  out  the 
magnitude,  yielding  the  density  of  the  phase  of  each  FFT  bin. 

(4)  Using  numerical  Integration,  obtain  the  mean  and  variance  of  the 
bearing  estimate  from  the  phase  density. 

The  density  of  the  phase  in  the  k1^  FFT  bln,  was  shown  to  be 
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1  ♦  erf (0cos 


(D“3> 


D-2 


erffx) 


dx,  x  2  0 


erf(x)  ■  -erf  (|*|),  x  <  0 


2  lEov^r 

3  "  Vir(W^) 


+  -  arg 


Wk  •  FFT  of  weight  vector  In  k“  FFT  bin. 

This  result  Is  equally  valid  In  the  case  where  the  both  broadband  and  narrowband 
directional  components  are  present.  From  (2],  when  both  are  present,  the  mean 
and  variance  of  the  weights  In  the  kth  FFT  bin  In  steady  state  is: 
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Therefore,  for  this  case 


Id  these  expressions 

M  ■  adaptive  filter  feedback  coefficient 


power  In  narrowband  directional  component 


variance  of  broadband  directional  component 


variance  of  noise 


M  ■  no.  of  adaptive  filter  taps  -  no.  of  FFT  bins 


ICflc.WJI 


T9  ■  algorithm  sample  rate 

u>0  ■  radian  frequency  of  narrowband  component 

T1  -  i 8in  *a  ’  AT.  ’  iaUx*b  (D~S 

Tj  ■  delay  In  narrowband  signal  between  array  phase  centers 
XT  ■  delay  in  broadband  signal  between  array  phase  centers 

O 

A  ■  Integer 

d  “  distance  between  array  phase  centers 
c  -  speed  of  sound 

#s  -  bearing  angle  of  narrowband  source 
■  bearing  angle  of  broadband  source 

It  is  assumed  that  the  directional  sources  are  In  the  same  horizontal  plane  as  the 
array  phase  centers. 

The  bearing  estimate  11]  Is  then 


i  ■  sin 
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where  a  Is  the  phase  of  the  bln  with  the  largest  magnitude.  Therefore, 
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r  ^  i  r^2  2P''v" 

Var  I  9  -  E  9  *  -  E  9 
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The  density,  f#  ,  is  that  of  the  phase  In  the  bln  with  the  largest  magnitude. 
While  these  lnt^$ftfis  cannot  be  evaluated  explicitly,  they  can  be  evaluated  numer¬ 
ically  on  the  computer.  Care  must  be  take?  In  the  Interpretation  of  these  results 
due  to  the  modulo  2*  assumption  Inherent  In  the  derivation  of  the  phase  density 
(see  pages  15  and  22  of  Appendix  A). 

The  density  given  in  (0-3)  Is  even  about  (•  -  y )  -  0,  so  that 
E(e  -V)  -  E[#]->  -  0 


E [*]  -V 


Since  the  estimate  of  Interest  la  the  bearing  of  the  narrowband  source,  the  "correct 
phase  angle  Is  eg  -  .  Expression  (20)  can  be  written  as 


a  rctan 
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with  X  the  signal- to- Interfere  nee  ratio  In  the  FFT  bln,  given  by 
I  G(k,  ->_)  1 2  <t\ 

- - Z  •. 

(M-A)  <r‘ 
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The  bias  In  the  phase  angle  Is  E(S  -  Tj)  «■  v  -  u>0  Tj.  By  a  trigonometric  identity 
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tan  r.  -  tan  V 


♦  tan  tan  > 


Substituting  (24)  and  (22)  and  solving  for  u>n  r 
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This  gives  the  bias  In  the  phase  angle  of  the  FFT  bin  containing  the  signal  as  a 
(Unction  of  the  time  delays  associated  with  the  broadband  and  narrowband  signals. 
Figure  D-l  shows  this  (Unction  plotted  against  the  difference  (ugT.  -  2*  to  MAk) 
for  various  values  of  the  signal-to-lnterference  ratio,  X  .  For  angles  from  -0°  to 
-180°,  it  can  be  seen  from  (25)  that  the  (Unction  is  odd  about  0°. 

When  the  apparent  bearing  angle  of  the  target  is  small,  then  the  bearing  esti¬ 
mate  is 
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Then  the  bias  in  the  narrowband  bearing  estimate  is 

B  -  E  j"d  -6  «  -tSt  >  -  9 
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The  first  term  is  a  bias  occurring  when  the  narrowband  component  is  not  bin  cen¬ 
tered,  and  becomes  negligible  as  the  resolution  of  the  FFT  increases.  The  second 
term  is  just  expression  p- 17)  multiplied  by  the  factor  (MTsc/2*kd). 

To  obtain  some  feel  for  the  magnitude  of  the  bias  to  be  expected  In  the  narrow- 
band  bearing  estimate  due  to  broadband  interference,  consider  an  array  with  75 
feet  between  phase  centers,  operating  on  a  600  Hz  sinusoid,  which  is  bln  centered 
In  the  FFT.  In  this  case, 


/MTs  c' 


0.0178 


t 


t{uit  is,  a  10°  bias  in  phase  a,  represents  only  a  0.  178°  bias  in  bearing.  Then  to 
interpret  Figure  D-l  for  this  array,  the  scales  on  the  axes  are  multlpled  by 
0.0178.  Hence,  for  an  interference  3dB  below  the  narrowband  component  in  the 
FFT  bin,  the  bias  is  at  most  (0.0178  x  30°)  -  0.531°,  while  for  an  interference 
12  dB  down,  the  bias  is  at  most  0.0637°. 
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Figure  D-l.  Bio*  in  Phase  Estimate  lor  Adaptive  Tracker  with  Broadband  and 
Narrowband  Targets 
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The  following  observations  should  be  noted: 

1.  The  apparent  Immunity  of  the  estimate  to  large  bearing  bias  errors  Is 
primarily  due  to  the  fact  that  the  bearing  estimate  la  ambiguous  be¬ 
yond  x*e/u>0d  radians,  so  all  estimates  are  Interpreted  In  this 
range. 

2.  When  the  signal  to  Interference  ratio  In  dB  la  negative,  the  bias  is 

between  (1/2)  (  "  0  ^  and  (0^  -  9^.  For  even  moderately  positive 

signal  to  Interference  ratios,  the  bias  la  considerably  smaller  than 
(1/2)  <e8  -  9^. 

3.  This  analysis  has  not  Included  any  array  response.  To  Include  array 
response,  the  relative  array  gains  for  the  narrowband  and  broadband 
components  would  be  Incorporated  In  the  signal  to  interference  ratio, 

X. 

4.  There  la  no  reason  to  expect  the  behavior  of  this  narrowband  estimator 
to  differ  significantly  from  any  other  narrowband  split  beam  estimator 
that  extracts  a  bearing  estimate  from  the  phase  difference  In  an  FFT 
bln  containing  the  signal. 

5.  From  'D-19),  It  can  be  seen  that  a  bias  occurs  when  ever  the  line 
frequency  Is  not  bln  centered.  This  bias  has  two  components,  one 
due  to  using  the  bln  center  as  the  frequency  estimate,  and  one  due  to 
a  bias  of  the  phase  of  the  FFT  bin.  This  bias  exists  even  when  the 
broadband  and  narrowband  sources  are  colocated,  as  In  the  case  of  a 
targe  t  with  the  broadband  and  narrowband  radiated  noise.  Increasing 
the  resolution  of  the  FFT  reduces  this  bias. 

Figures  D-2  through  D-4  show  the  variance  of  the  bearing  estimate  for  the  nar¬ 
rowband  adaptive  tracker  with  a  target  having  both  broadband  and  narrowband 
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Fipirr  D-2.  Standard  (Vrution  of  Bearing  Error  Vcnui  Narrowband  SNR 
for  Various  Broadband  SNR 
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component  a,  determined  by  the  numerical  method  above.  The  narrowband  compo¬ 
nent  la  bin  centered,  ao  that  the  eatlmate  la  unbiaaed.  The  almulatloo  utilizes  a  0 
to  256  Ha  band  with  a  2  Hz  FFT  resolution,  and  an  array  with  150  ft  between  phaae 
centers.  The  variance  la  plotted  aa  a  function  of  the  ratio  of  narrowband  algnal 
power  to  background  nolae  (narrowband  SNR)  for  varloua  value  a  of  broadband  alg¬ 
nal  power  to  nolae  ratio  (broadband  SNR).  The  plot  a  only  extend  over  the  range 
of  parametera  for  which  the  filter  la  atable.  Plata  for  p  -  2'12,  2~14,  and  2~ 16 
are  given  In  Figures  D-2,  D-3,  and  D-4,  respectively.  In  all  caaea,  the  target 
la  at  0°  azimuth.  Figure  2  shows  several  points  from  computer  simulations  at 
M  •  2* which  verify  the  analytical  results. 

Aa  noted  earlier,  the  analyses  leading  to  the  expressions  (D-4)  and  (D-5)  as¬ 
sume  that  the  correct  FFT  bin,  1.  e. ,  the  one  containing  the  narrowband  component, 
la  chosen  for  bearing  extraction.  This  selection  la  baaed  upon  the  FFT  bin  having 
the  largest  magnitude.  When  both  narrowband  and  broadband  components  are 
present,  the  total  broadband  energy  (algnal  plus  noise)  may  result  In  choice  of  the 
Incorrect  bln.  While  the  narrowband  tracker  does  show  the  ability  to  determine 
bearing  from  broadband  signals,  this  Incorrect  selection  does  change  the  statistics 
of  the  estimate,  as  the  estimate  depends  on  bln  number,  via  (D-10).  In  the  simula¬ 
tion  results  shown  In  Figure  D-3,  the  algorithm  baa  been  forced  to  choose  the  cor¬ 
rect  bln,  regardless  of  narrowband  to  broadband  power  ratio,  in  order  to  verify 
the  results  of  the  analysis.  Therefore,  when  the  narrowband  to  total  broadband 
power  ratio  in  the  FFT  bln  is  less  than  about  10  dfi,  the  curves  are  not  valid  un¬ 
less  some  other  means  of  determining  the  correct  bin  is  available.  For  example, 
an  additional  FFT  with  post- FFT  Integration  might  be  uaed  to  Identify  the  bin  con¬ 
taining  the  narrowband  line. 
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The  analysis  described  above  can  also  be  used  to  assess  the  performance  of  the 
narrowband  adaptive  tracker  with  a  narrowband  target  and  a  broadband  Interference 
in  a  different  direction.  Figures  D-5  through  D-8  show  the  effect  of  a  broadband 
Interference  at  2°  and  4°  in  bearing,  on  the  bearing  estimate  for  a  narrowband  tar¬ 
get  at  0° .  The  configuration  of  the  tracker  Is  Identical  to  that  used  In  Figure  D-3 
with  m-  2"14.  Figures  D-5  and  D-6  show  the  standard  deviation  and  bias  of  the 
estimate  for  the  interference  at  2°,  while  Figures  D-7  and  D-8  give  the  same  In¬ 
formation  for  an  Interference  at  4°.  As  would  be  expected  1 21,  when  the  inter¬ 
ference  power  In  the  FFT  bln  Is  significant  relative  to  the  power  of  the  target  line, 
the  estimate  la  biased  toward  the  Interference.  Additionally,  Figures  D-5  and  D-7 
show  that  for  sufficiently  high  SNR,  the  presence  of  a  broadband  interference  may 
Increase  the  variance  of  the  estimate.  This  Increase  In  variance  appears  to  occur 
in  a  range  of  SNR  where  bias  due  to  the  Interference  is  no  longer  significant.  At 
lower  SNR,  the  variance  of  the  estimate  is  reduced  with  Increasing  Interference 
power,  while  the  bias  becomes  significant  (l.e. ,  the  output  is  a  "good"  estimate 
of  the  lnterfsrence  bearing).  The  SNR  at  which  the  transition  from  varianoe  de¬ 
creasing  with  LVR  to  variance  increasing  with  INR  appears  to  decrease  with  In¬ 
creasing  separation  between  target  and  Interference. 

Figures  D-9  and  D-ll  show  the  variance  of  the  estimate  for  the  same  tracker 
configuration  when  the  target  Is  at  2°  and  4°  respectively,  and  the  Interference  Is 
at  0*  and  Figures  D- 10  and  D- 12  show  the  bias  for  the  same  situation.  These  curves 
show  the  results  of  simulations  of  the  device  under  the  same  Input  conditions. 
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Figure  D-7.  Effect  of  BB  Interference  on  MB  Tracker  with  NB  Target  (Std  Deviation  of 
Bearing  Error  va  SNR  for  Various  INR) 
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Figure  D-10.  Effect  of  BB  Interference  on  NB  Tracker  with  NB  Target  (Bias  in  Bearing 
Estimate  n  SNR  for  Various  INR) 
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Fipue  D-ll.  Effect  of  BB  [nterfcrence  on  NB  Tracker  (NB  Target)  (Std  Deviation  of 
Bearing  Error  vi  SNR  for  Vannus  I  NR) 
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Appendix  E 

Comparison  of  an  Adytive  Tracker  and  a  Split  Beam  Tracker 
with  Incorrect  a  Priori  Assumptions  on  Input  Spectra 

The  broadband  aduptive  tracker^1®  has  been  suggested  as  an  alternative 
to  conventional  split  beam  correlator  trackers  because  of  several  potential 
advantages  stemming  from  the  adaptive  characteristics.  One  of  these 
possible  advantages  is  that  the  adaptive  filter  does  not  require  a  priori 
knowledge  of  the  signal  and  noko  statistics,  but  estimates  these  dynamically 
during  operation.  Conventional  split  beam  trackers,  on  the  other  hand, 
require  the  input  signal  and  noise  spectra  in  order  to  optimize  their  per¬ 
formance  i see,  for  example,  161  or  Ul).  The  goal  of  this  analysis  is  to 
determine,  at  least  for  some  simple  cases,  whether  or  not  the  adaptive 
tracker  (ADT>  is  more  tolerant  of  sp**ctral  variation  about  those  assumed 

a  priori  than  the  split  beam  tracker  <  SB T )_, - 

The  approach  is  as  follows.  A  signal  and  noise  srpectrum  is  assumed 
a  priori .  and  the  weighting  Alter  for  the  SBT  is  set  using  these  spectra. 
Then,  with  the  input  spectra  matching  the  a  priori  assumptions  the 
variance  of  the  bearing  estimate  for  the  ADT  and  SBT  .re  computed 
analytically.  Either  the  input  noise  or  signal  spectrum  is  then  varied 
without  changing  the  weighting  filter  of  the  SBT.  and  the  variances 
recomputed.  The  ratio  of  the  variances  with  correct  and  incorrect  a  priori 
kw-’-dge  gives  the  change  in  performance  of  the  two  trackers.  These 


a 


«d  to  determine  which  of  the  two  is  less  susceptible  to 


changing  inputs.  In  this  report,  the  assumed  signal  and  noise  spectra  are 
flat  over  the  input  band  of  the  tracker.  The  case  in  which  the  actual  spectra 
tire  of  some  general  shape  is  so  complicated  that  it  yields  little  insight. 

The  most  detailed  results  are  for  the  specific  examples  in  which  actual 

spectra  differ  from  the  assumed  only  in  level  in  some  portion  of  the  input  band. 

Analysis: 

Let  the  split  inputs  to  the  tracker  be 


x.(t)  -  S(t  -  —  sin  0  )  ♦  nT  (t) 
L  C  L 


xR(t)  =  S(t)  ♦  nR(t) 


f 


(E-l) 


where  S(t),  n^(t),  and  nR(t)  are  independent,  zero  mean  gauaaian  random  proc¬ 
esses  with  power  spectral  densities  SB(~)  and  Sq(w)  (q^  and  nR  have  identical 


power  spectra) .  Also-,  let 


d  -  distance  between  half  array  phase  centers 
c  *  speed  of  sound 

0  =  angle  between  plane  wavefront  and  axis  of  phase  centers. 


In  Appendix  4  of  (1] ,  the  mean  and  variance  of  the  transfer  function  of 
a  continuous  adaptive  Alter  is  derived  under  the  condition  of  equal  signal 
and  noise  bands.  This  was  done  by  resolving  the  inputs  into  discrete  fre¬ 
quency  bins,  and  applying  a  single  tap  complex  adaptive  filter  in  each  bin. 
as  discussed  in  Appendix  3  of  ( 1 ) .  The  transfer  functions  are  obtained  as 
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the  limit  of  the  discrete  case.  This  approach  depends  upon  the  independence 
of  the  individual  frequency  bins.  The  resulting  transfer  function  an  be 
used  to  obtain  the  impulse  response  of  the  adaptive  filter,  h(  T)* 

The  angle,  e  .  is  to  be  estimated  by  determining  the  value  of  r  for 
which  the  impulse  response.  h(r),  is  a  maximum.  This  is  an  estimate  of  the 
delay  between  arrivals  of  the  signal  wavefront  at  the  half  array  phase  cen¬ 
ters.  and  can  be  converted  to  a  bearing  estimate.  d  ,  using 

t=  -  sin  d  .  (E-2) 

c 


The  value  of  0  for  which  h(  t  )  is  peaked  corresponds  to  the  value  for 
which  the  derivative  of  h(  r )  with  respect  to  6  is  zero.  In  the  neighborhood 
of  the  zero  crossing  the  derivative  of  the  Impulse  response  is  approximately 
linear.  The  fluctuations  in  the  derivrtive  of  the  transfer  function  then  map 

- —  through  the  linear  function  to  provide  an  estimate  of  the  mean  square  error 

in  the  angle  estimate.  Thus  the  mean  value  of  the  derivative  of  the  trans¬ 
fer  function,  the  standard  deviation  of  the  derivative  of  the  transfer  func¬ 
tion.  and  the  slope  of  the  mean  of  the  derivative  of  the  transfer  function 
are  required,  at  the  point  where  0  =  8.  Then  the  errors  in  the  transfer 
function  can  be  mapped  into  the  errors  in  one's  aoility  to  extract  the  peak 
of  the  transfer  function  as  follows: 


Var 
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This  approach  Is  motivated  by  the  treatment  of  BDI  split  beam  trackers  by 

MacDonald  and  Schultheiss  in  Reference  [  2] . 

We  want  to  compare  the  adaptive  tracker  with  the  conventional  split 

beam  tracker,  which  can  be  done  using  MacDonald  and  Schultheiss’ 

16] 

results.  These  results  require  a  priori  assumptions  regardiig  the 
signal  and  noise  spectra.  The  a  priori  signal  and  noise  spectra  will  be 

denoted  and  Sna(w),  respectively,  not  necessarily  equal  to  the 

actual  spectra. 

Hero,  we  will  examine  the  effect  of  incorrect  a  priori  assumptions  on 
the  performance  of  the  split  beam  tracker  in  several  cases.  Thi3  will  be 
done  by  comparing  the  variance  of  the  bearing  estimate*  under  actual 
conditions  with  the  case  when  a  priori  assumptions  are  true.  The  per¬ 
formance  of  the  adaptive  filter  under  both  conditions  will  aU>  be  computed. 
In  order  to  maintain  the  same  time  constant  for  the  ADT .  the  power  at  the 
input  is  normalized  (equivalent  to  an  AGC)  such  that 

f  "  1  Ss(w)  *  Sn(w)1  a  l  s,a(w)  ♦  *na(u,)1  **  (E-4) 
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while  the  actual  spectra  are  Ln(u»)  =  Ln#(w)  and 
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These  spectra  are  illustrated  in  Figure  E-l,  where  Kj  is  selected  to  satisfy 
l£-4)t  so  that 
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The  results  cited  in  reference  [1]  lead  to  the  following  statistics  for  the 
transfer  function  of  the  adaptive  filter. 


This  is  to  be  compared  with  the  split  beam  tracker  under  the  same  cir¬ 


cumstances.  From  MacDonald  and  Schutheiss 


beam  tracker  (SBT)  bearing  estimate  is  given  by 


where  H(^)  is  a  weighting  function  given  by 


so  factoring  this  out  of  (E-22)  yields 
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At  low  SNR  we  get 


In  the  special  case  of  Lgl  =  0, 


which  gives,  for  small  SNR. 


Comparing  (E-15)  with  (E-25),  and  (E-17)  with  (E-27)  for  the  case 

when  L  ,  3  0.  it  can  be  seen  that  at  low  SNR,  the  SBT  and  the  adaptive 
si 

tracker  degrade  in  the  same  way  from  the  case  when  the  a  priori  assumptions 
are  true. 


Case  2:  Noise  spectrum  differing  from  a  priori  assumptions. 

Here,  the  a  priori  assumptions  are  the  same  as  in  case  1,  but  the 
actual  spectra  are.  as  shown  in  Figure  E-l, 
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0.  otherwise 


This  is  a  straightforward  extension  of  case  1.  giving 
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Were  the  a  priori  assumptions  valid,  the  variance  would  be  Vq.  given  in 
tE-13),  above,  so 
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As  in  Case  1.  we  are  interested  mainly  in  the  low  SNR  case  where 

L  «L  and  L  «L  , ,  so 
s  n  s  n  l 
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In  this  case.  K,  is  selected  such  that 
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For  the  SBT.  we  obtain 
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From  (E-37),  it  c.  n  be  seen  that  for  all  values  of  K^,  leas  than  unity,  the 
variance  of  the  SBT  bearing  estimate  la  monotonically  Increasing  In  Kn>  Hence 
when  the  actual  noise  in  the  band  w  ^  to  &  is  less  than  assumed,  the  SBT  performs 
better,  while  if  there  is  more  noise.  It  degrades.  It  should  also  be  observed  that 
incorrect  a  priori  assumptions  as  to  the  noise  spectrum  alfects  only  the  variance 
of  the  correlation  tunction,  not  the  mean,  while  affecting  both  the  mean  and  the 
variance  of  the  ADT  weights.  Figure  E-2  shows  the  performance  factor,  K^, 
plotted  as  a  function  of  for  various  values  of  K  .  When  is  greater  than 
unity  the  variance  increases,  while  It  decreases  for  less  than  one. 

The  parametric  behavior  of  the  AOT  Is  not  so  obvious  from  (E-36).  By 

definition,  the  variance  of  the  estimate  is  V  when  K  *  1,  and  it  can  be  shown  that 

on 
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This  is  positive  for  all  values  of  K  t  (0.1],  so.  at  least  in  the  vicinity  of 

Kn  ■  1  (a  priori,  assumptions  approximately  true),  the  variance  of  tho  ADT 

estimate  increases  with  K  .  Further 
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E-15 


F *urr  E-2.  Chingf  in  Variance  of  Split  Beam*  Tracker  Bearing  Estimate  with  Incorrect 
A  priori  None  Spectrum  Input 
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so  for  large  K^.  the  ADT  degrades  from  its  performance  when  a  priori 
assumptions  are  true.  Finally,  as  Kn— ►  0, 

Var(0)  —►  0 


Figure  E-3  show  the  performance  factor  for  the  ADT 
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plotted  as  a  function  of  K  for  various  values  of  K  .  It  can  be  seen  that, 

n  j 

as  in  the  case  of  the  SBT.  the  variance  of  the  estimate  increases  with  K  . 

n 

In  order  to  compare  the  way  in  which  the  two  trackers  degrade,  con¬ 
sider  the  ratio  of  their  variances  at  low  SNR. 
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Figure  E-J.  Change  in  Variance  of  Adaptive  Tracker  Bearing  Estimator  with  Changes 
in  NoiK  Spectrum. 


For  K  <  1,  when  both  trackers  degrade  trom  the  variance  with  correct  a  priori 
n 

assumptions,  D<  1  indicates  that  the  ADT  is  less  susceptible  to  degradation. 

When  >  1,  where  both  traders  improve  performance,  D<  1  indicates  that  the 
ADT  has  improved  more  than  the  SBT.  It  can  be  shown  that 

(1)  If  Kn  S  1.  D  <  1 

»  *».  D  —  0 

-  0.  D  -  0 

>  0 

n=  1 

Therefore,  for  Kn  <  1  and  for  large  K^.  the  adaptive  tracker  degrades  less 

severely  than  the  SBT.  but  by  iE-4),  the  ADT  is  stgwrlor  at  least  for  some 

values  of  K  greater  than,  but  close  to.  one. 
n 

Figure  E-4  shows  the  performance  ratio,  D,  plotted  as  a  (unction  of  for 

various  values  of  K  .  The  value  of  D  is  less  than  one  except  in  a  region 

extending  from  K  -  1  to  some  value  of  K  ■  K  .  where  K  >1.  The 
n  n  no  n 

value  of  Kno  appears  to  be  Inversely  proportional  to  K  .  In  the  region 

where  the  SBT  is  superior.  It's  advantage  is  slight,  while  the  ADT  may  be 

significantly  better  than  the  SBT  for  values  of  K  <1  and  K  >  K  .  * 

n  n  no 

Figure  E-5  shows  the  magnitude  and  location  of  the  maximum  of  D  as  a  func¬ 
tion  of  Ku .  At  most,  the  ADT  bearing  estimate  is  degraded  0.53  dB  from 
that  of  the  SBT.  while  it  is  clear  frou  Figure  E-5  that  the  ADT  offer 
very  significant  advantages  for  other  values  of  K  .  For  reasonable 
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values  of  K  .the  ACT  Is  superior  except  when  the  noise  power  in  the  band 


is  increased  from  0  to  3  dB,  in  which  case,  the  SBT  may  provide 

n  a 

slightly  better  results. 

It  is  possible  to  develop  expressions  for  the  relative  degradation  of  the  two 
trackers  when  the  actual  signal  or  noise  spectrum  is  arbitrary,  but  assumed 
white  a  priori.  These  derivations  do  not,  however,  provide  much  insight  but  are 
included  below. 

Conclusions 

Two  simple  cases  of  incorrect  a  priori  assumptions  on  input  spectra 
and  their  effect  upon  both  the  SBT  nnd  the  adaptive  tracker  have  been 
considered  at  low  signal  to  noise  ratio.  In  the  case  where  signal  assump¬ 
tions  were  incorrect,  the  ADT  and  SBT  degraded  identically,  primarily 
due  to  the  low  SNR.  In  the  case  where  the  noise  spectrum  differed  from 
that  assumed,  the  ADT  offered  significant  immunity  to  degradation  relative 
to  the  SBT  for  most  spectral  variations  from  that  assumed.  For  the  limited 
cases  where  the  SBT  degrades  less  t.ian  the  ADT.  its  advantage  is  at  moat 
0.5  dB.  The  ADT  appears  to  be  able  to  better  take  advantage  of  reduced 
noise  in  some  portion  of  the  input  band,  and  to  reduce  noise  levels  when  a 
large  Increase  occurs  in  part  of  the  band. 
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STATISTICS  OF  FREQUENCY  DOMAIN  WEIGHT  VECTOR  FOR  NON-FIAT 
INPUT  SKjNAL  AND  NOISE 

In  Appendix  IV  of  reference  [1].  the  statistics  of  a  frequency  domain 
weight  vector  of  a  continuous  adaptive  filter  were  developed  for  signal  and 
noise  spectra  that  are  flat  between  two  frequencies,  wT  and  u  .  that  is 
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The  statistics  of  the  transfer  function  of  the  adaptive  filter.  H(ui).  in 
these  cases  are 
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In  this  appendix,  these  statistics  are  developed  for  signal  spectrum. 

S  (w),  and  noise  spectrum  S  (w).  symmetric  and  band  limited  to  fw.  .  w  )  as 

9  n  L  U 

follows.  A  stepwise  approximation  to  the  signal  spectrum  is  constructed  as 
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Using  the  results  from  reference  111  on  each  of  these  intervals,  the 
statistics  of  the  adaptive  filter  transfer  function  are 
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To  compute  the  variance  of  the  bearing  error  for  the  adaptive  tracker,  we 
require 
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Again  applying  the  results  of  Appendix  IV  of  reference  (  l  ] ,  we  get 
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Note  that  these  are  both  Riemann  sums,  with 
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COMPARISON  OF  THE  ADT  AND  SBT  WITH  ASSUMED  WHITE  SIGNAL 
AND  NOISE,  AND  GENERAL  SIGNAL  SPECTRUM 


If  the  prior  assumptions  were  true,  tne  variance  would  be 
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so  factoring  this  out  of  the  above  gives 
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We  can  use  equations  (IS)  and  (19)  to  obtain  tne  variance  for  the  SBT 


under  the  same  conditions,  where 
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Then,  at  low  SNR 


At  low  SNR.  then,  the  adaptive  tracker  and  SBT  degrade  the  same  way 
under  incorrect  a  priori  signal  assumptions. 


COMPARISON  OF  THE  ADT  AND  SBT  WITH  ASSUMED  WHITE  SIGNAL 
AND  NOBE  AND  GENERAL  NOISE  SPECTRUM 

Let  the  a  priori  assumptions  be  as  in  Case  II.  but  let 
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and  let  Sn(  ui)  be  a  symmetric  spectrum  band  limited  to  I  «*>  I  s  ‘‘'a  • 


Proceeding  .is  in  Appendix  II.  under  small  SNR. 
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Similarly,  the  variance  for  the  SBT  is  as  follows: 
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Figure  CM2.  Effect  of  BB  Interference  on  NB  Tracker  (NB  Target)  (Biaa  in  Bearing 
Estimate  vs  SNR  for  Various  INR) 


Broadband  Bearing  Estimating  in  the  Presence  of  an  Interference 


During  the  first  phase  of  this  study,  the  performance  of  the  adaptive 
tracker  with  a  broadband  input  signal  was  analyzed.  In  this  section,  the 
ground  work  for  an  analysis  of  the  effect  of  a  single  broadband  interference 
plus  ambient  noise  on  the  bearing  estimate  for  a  broadband  target  is  developed. 
Here  a  broadband  interference  refers  to  a  directional  point  source  radiating 
broadband  noise,  generally  separated  from  the  target  in  azimuth.  The  ap¬ 
proach  to  the  analysis  of  the  adaptive  tracker  performance  in  the  presence 
of  an  interference  will  be  a  direct  extension  of  the  method  used  in  the  analy¬ 
sis  for  broadband  interference  in  the  first  phase  of  this  study.  The  analysis 
will  initially  utilize  a  continuous  version  of  the  adaptive  filter,  as  shown  in 
Figure  F-l.  By  analogy  with  the  frequency  domain  model  of  Figure  A-2,  it 
is  possible  to  determine  the  mean  and  variance  of  the  continuous  weight 


function  in  the  presence  of  interference.  As  before,  the  time  delay  estimate 
between  split  arrays  is  extracted  by  locating  the  peak  of  the  weight  function. 

A  crucial  difference  between  the  signal  only  and  signal  plus  interference 
case  is  that  in  the  latter,  the  weight  function  may  have  two  peaks,  one  at 
the  signal  delay  and  the  other  corresponding  to  the  delay  in  the  interfer¬ 
ence  between  split  arrays.  Moreover,  there  may  be  only  a  single  peak 
located  between  the  two  delays.  When  two  peaks  occur,  the  estimate  will  be 
unbiased  if  we  select  the  correct  (signal)  peak.  If  only  a  single  peak  occurs 
between  the  two  delays,  then  the  estimate  will  be  biased. 


As  in  Appendix  IV  of  reference  [1]  .  and  in  [6],  the  variance  in  the  deter¬ 


mination  of  the  location  of  the  peak  is 
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where  h(t)  is  the  weight  function,  t  is  the  estimate  of  the  peak  location,  and 
Tp  is  the  location  of  the  peak  in  the  mean  weight  function.  This  result  can 
he  easily  converted  to  a  bearing  error  for  a  given  array. 

Just  as  in  the  analysis  of  the  signal  only  broadband  case,  the  array  is 
treated  as  two  omnidirectional  point  sensors  located  at  the  phase  centers  of 
the  split  arrays.  This  model  is  valid  when  the  signal  and  interference  are 
near  the  MR  A.  so  that  the  response  of  the  half  arrays  are  nearly  invariant 
with  frequency. 


Determination  of  Mean  and  Variance  of  Weight  Vector  in  the  Presence  of 
Broadband  Interference 

For  the  broad>and  signal  plus  interference,  the  time  domain  inputs  to 
the  adaptive  tracker  can  be  modeled  as 
x(t)  =  s(t)  ♦  n  (t) 

(F-2) 

d(t)  *  s(t  -  ATs)  ♦  i(t  -yTj)  ♦  nd(t) 
where  s(t),  i(t),  nT(t).  and  n^it)  are  zero  mean,  stationary  random  pro¬ 
cesses.  uncorrelated  in  time  and  independent  of  each  other.  Let  the  delay 

between  the  two  inputs  for  signal  and  interference.  A T  andyT  .  be  multiples 

s  s 

of  the  sampling  interval.  T#. 

The  mean  and  variance  of  the  weight  vector  can  be  analyzed  using  the 
frequency  domain  model  for  the  adaptive  filter  shown  in  Figure  F-l.  An  FFT 
is  performed  on  both  x(t)  and  d(t)  every  RT9  seconds,  so  that  each  frequency 
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domain  filter  iterates  every  RT  seconds.  With  M  the  FFT  size.  R  =  M  oorre- 
s ponds  to  FFT  processing  without  overlapping  or  gaps  between  time  windows. 
R  <  M  corresponds  to  overlapping  FFTs  and  R  >  M  Indicates  gap  processing. 
The  ith  input  to  the  adaptive  filter  in  the  kth  FFT  bln  is  given  by 
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Lsing  the  results  of  reference  (1]  ,  Appendix  III,  the  mean  square  value  of 
the  weight  in  the  kth  FFT  bin  is  defined  by 
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The  single  sum  is  treated  first . 


We  also  need  |E[Dk(m)  X  (m)l  Consider  the  signal  contribution. 
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By  analogy  with  the  signal  term  , 


E[|Dk(m)Xk(m)|2]  =  M2(js2  ♦  Oj2  ♦  z^) 


4  2t  4  2-t  .  2 

2  'll  ^  2  ll  yk 

♦  (M- A)  c  e  M  ♦  (M-y )  o,  e 


J  XT  Ak  2 

♦  S(M.A)  e  a.  *  S(M.y)  e 


, 2»  .  ! 0 
J  TT  *  2  l  2  (F-101 

a.  1 


We  also  need 
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The  double  sum  term  in  E[W(n«-l)]  .  involves  E[F(m)F*(q)J  which  is 
treated  as  follows: 
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In  the  double  sum.  q  <  m ,  so  that 
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Thus,  since  both  p  and  v  ranfce  over  indices  that  are  greater  than  m  and  q, 


E[Fkfm)F*(q)]  =  E  TT  TT  IK  Xk(p)  2](1  -V  Xk(q)  2] 


I 


p=m*  1  v»m*l 


1 


(F-18 


m 


•  E  Dk(m)Xk*(ro)Dk*(q)Xk(q)  1 1“  »* ;  xk<v>  21 

'  v*q*l 


The  flrst  expectation  was  evaluated  previously  in  equation  (F-fl).  The  second 
expectation.  since  q  is  always  leas  than  m.  is 


(F-24) 


Taking  the  limit  ee  n  *  •  yields  the  steady  state  mean  square  weight. 
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Then  the  steady  etete  weight  variance  is 
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This  agree*  with  the  result  for  the  cese  of  Uj  «  0  end  &*0  derived  in 
Reference  [1],  Appendix  III. 

We  are  interested  in  the  performance  of  the  time  domain  adaptive  tracker. 


for  which  windowing  effects  do  not  occur.  Then  in  this  case. 
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This  is  the  same  form  as  the  variance  without  the  interference  but  with  the 
signal  power  replaced  by  the  sum  of  the  signal  and  interference  powers. 


Determination  of  Bias  and  Bearing  Estimate  Variance  for  Broadband 
Tracker  with  broadband  Interference 

As  in  reference  [  1  ] ,  it  is  assumed  that  the  mean  bearing  estimate  occurs 
at  the  peak  of  the  mean  weight  function.  E(h(r)  1 .  or  equivalently,  at  the 
zero  crossing  of  the  derivative.  E[3h(r)/drl .  corresponding  to  this  peak. 
Suppose  we  assume  that  the  input  signals  and  noise  are  flat  for  |w|c[u>Lwu) 
and  xero  outside  that  range.  Then  using  expression  (F-28),  the  mean  trans¬ 
fer  function  of  the  continuous  adaptive  Alter  can  be  seen  to  be 
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The  derivative  of  the  mean  weight  function  in  the  time  domain  is  then 
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We  are  interested  in  Anding  the  location  of  the  xero  crossing  of  this  function 
corresponding  to  the  maximum  of  the  weight  function.  Although  this  cannot 


be  determined  explicitly,  it  can  be  found  numerically  on  the  computer.  For 
convenience,  we  will  consider  the  lowpaas  case,  where  =  0,  and 
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where 


t  3  -sin?  and  i.  *  -sin?, 
s  c  s  I  c  I 


To  compute  the  variance  of  the  estimate  from  (F-l),  we  also  need 
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Finally,  we  need 
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since  we  will  choose  r  such  that 


Note  this  expression  is  actually  valid  only  for  values  of  tp  such  that 


The  actual  variance  is  obtained  by  substituting  the  value  of  rp  from  the 

numerical  maximization  of  (F-31)  in  this  expression. 

Bearing  Estimation  with  a  Discrete  Adaptive  Tracker  in  the  Presence  of 
an  Interference 

The  analysis  in  this  section  has  followed  the  same  course  as  that  of 
Appendix  IV  of  ID  .  using  a  continuous  adaptive  filter.  As  pointed  out  in 
Appendix  V  of  that  reference,  however,  in  practice  the  adaptive  filter  is  dis¬ 
crete  in  time  and  of  finite  length,  and  the  peak  of  the  impulse  response  must 
be  found  by  interpolating  between  discrete  sample  points.  This  interpola¬ 
tion  increases  the  variance  of  the  bearing  estimate  with  respect  to  that  of 
the  continuous  case.  The  interpolation  was  of  the  type  in  which  an  inter¬ 
polated  impulse  response  was  computed  as 
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hI(t)  *  2  h(mTs)  f(t  -  mTs)  (F-38) 

m=0 


where  h(mTs)  are  the  discrete  filter  weights  and  f(t)  is  a  interpolation 
function .  such  as 


f(t)  * 

2  Bt 


The  algorithm  then  determines  the  peak  of  hj(t)  as  the  estimate  of  time 
delay  between  half  array  phase  centers. 


F- 19 


Now.  assume  that  the  inputs  to  the  adaptive  filter  are  bandlimited  to  some 


bandwidth.  B^.,  so  that 

H(k)  =  0  for  k  >  (J-l)  =  lBsTg J  J  (F-42) 

where  l X J  is  the  largest  integer  less  than  or  equal  to  X.  Further,  in  order 
for  h(mTs)  to  be  real,  we  require 

H(k)  =  H  *(M-k)  (F-43, 

Using  ( F- 4 1  > .  it  oar  readily  be  verified  that 

M-l  M-l  M-l 

Var1'2  [h’,(t>]  =  ^  Var  (H(k)l  ^  ^  ftt-m^) 

"  k=0  mj=0  m2=0 

i^*(m  -m,)k 

•r(t-m2Ts)em  *  ( F- 44) 

and 


We  want  (F-28)  and  (F-29)  to  satisfy  (F-42)  and  (F-43).  so 
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Then,  substituting  in  (F-44)  and  (F-45)  gives 
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and 


Then  the  variance  of  the  estimate,  given  t  .  is  found  by  substituting  (F-48) 
and  ( F-  49)  in  (F-39).  Given  any  interpolation  function.  f(t).  and  a  value 

for  t  .  the  Variance  of  the  delay  estimate  can  be  evaluated  numerically  on 

s 

the  computer.  The  value  of  t^.  the  location  of  the  peak  can  be  found 
numerically  from  (F-41). 


F- 


It  must  be  noted  here  that  this  result  applies  only  when  the  variance  of 


hj(t)  falls  in  the  linear  portion  of  h'j(t)  near  the  peak  of  E[hj(t)).  This  is 
equivalent  to  saying  that  the  largest  of  the  discrete  weights  must  be  in  the 
vicinity  of  the  true  peak.  This  will  be  true  when  the  mean  of  the  largest 
weight  is  large  in  comparison  to  the  variance  of  the  weights.  For  a  particu¬ 
lar  input  signal  to  noise  ratio,  this  can  be  assured  by  selecting  the  feedback 
coefficient  ,u.  sufficiently  small. 
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APPENDIX  G 


Appendix  G 

Tracking  of  Moving  Broadband  Target* 


Introduction 

In  (1).  the  variance  of  stationary  bearing  estimates,  using  the  adaptive 
filter  bearing  tracker,  was  compared  to  the  Cramer-Rao  lower  bound  on  the 
variance  of  any  unbiased  static  estimator.  The  comparison  showed  about  a 
0.5  dB  loss  for  the  adaptive  tracker.  As  a  further  comparison,  the  adaptive 
tracker  was  compared  to  the  split  beam  correlator  bearing  tracker.  In  (6] 
it  was  shown  that  the  ratio  of  the  variance  of  the  split  beam  correlator 
tracker  to  the  Cramer-Rao  lower  bound  is  given  by  [6  -  Eq  (33)]. 


where  M  =  total  number  of  hydrophones.  When  M  =  2.  as  when  comparing  to 
the  adaptive  tracker,  the  split  beam  correlator  achieves  the  lower  bound. 
Hence,  on  the  basis  of  static  estimation,  for  broadband  signals,  an  adaptive 
tracker  and  a  non-adaptive  split  beam  correlator  tracker  are  comparable. 
This  Appendix  addresses  operation  as  a  dynamic  device  (i.e. ,  tracking 
changing  bearing) .  It  Is  therefore  necessary  to  Investigate  the  dynamic 
tracking  capabilities  of  the  two  devices  and  make  a  comparison  of  each  in  a 
tracking  mode.  In  this  case  the  CR  bound  is  no  longer  applicable. 

The  variance  of  the  estimates  of  the  relative  time  delay  between  the  two 
half  arrays  is  computed  for  both  trackers  under  a  number  of  simplifying 
assumptions.  The  mean  weights  of  the  adaptive  filter  are  determined  using 
an  exponential  correlation  function  for  the  input  signal  correlation  for  a 


linearly  changing  delay.  The  variance  of  the  weights  with  signal  present  is 
approximated  by  the  noise  alone  variance.  These  two  expressions  are  used 
to  compute  the  variance  of  the  delay  estimate  in  the  absence  of  bias  errora. 
The  location  of  the  peak  of  the  weights  is  then  analytically  determined  and 
shown  to  lag  the  true  delay,  introducing  a  bias  into  the  estimate. 

For  purpoaes  of  comparison,  the  performance  of  a  first  order  discrete 
time  split  beam  bearing  deviation  indicator  tracker  is  Investigated.  The 
variance  of  this  tracker  is  determined  as  a  function  of  the  statistics  of  the 
correlator.  The  bias  error  is  determined  for  the  same  linearly  changing 
delay  as  for  the  adaptive  filter. 

The  two  devices  are  then  compared  on  the  basis  of  the  variance  of  their 
estimates  for  the  same  bias  error.  In  this  way.  each  device  is  designed  to 
achieve  the  same  dynamic  performance  with  some  resultant  fluctuating 
errors.  An  expression  is  then  derived  for  the  ratio  of  the  two  fluctuation 

error  variances. 

Adaptive  Filter  Bearing  Tracker  Model 

In  HI,  Appendix  VI,  an  expression  is  derived  for  the  mean  weights  of 
the  adaptive  tracker,  for  small  signal- to- noise  ratios  when  the  input  signal 
is  s  broadband  correlated  process  with  linearly  time-varying  delay,  buried 
in  uncorrelated  noise. 
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where 


adaptation  parameter 


„  2  2 
a  *  7n 


=  signal  and  noise  powers  respectively 


o  *  e  *  exponential  correlation  function  of  signal 

c  =  linear  rate  of  change  of  delay  between  two  phase 

centers  of  the  split  array 

It  is  shown  in  Addendum  G  I.  at  the  end  of  this  Appendix,  that,  for 
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tn  Addendum  GI.  st  the  end  of  this  Appendix,  it  Is  shown  thst  the  lag  In  the 
peak  of  the  weights  in  Eq  (G-3)  in  comparison  to  the  true  delay.  Is  given  by 
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Since  n  r  « 1  and  be  6  oc  1,  the  tap  lag  is  accurately  approximated  by 


(G-4) 


xc 


be  6  -  m  r 


-hrs 


be  6 


4 


If,  in  addition,  x*bc  6  ,  then 
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Eqs  (G-4)  and  (G-S)  have  been  shown  in  good  agreement  with  the  lag 
displayed  by  numerical  evaluation  of  Eq  (GI-4). 

It  is  shown  in  Addendum  Gil  that,  for  (1  -uon2)  a  pCA 
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Again,  using  the  condition  that  j  -  cl  *  0,  90  that 
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Eq  (G-7)  doer  not  peak  at  the  unbiased  value  of  delay,  j  =  c(n-l),  but  also 
displays  a  delay  in  the  location  of  the  peak .  Precise  expressions  for  the 
delay  will  not  be  investigated  for  this  case. 


Split  Beam  Tracker  Model 

Figure  o-l  shows  s  simplified  version  of  s  time  domain  split  beam  tracker 
for  two  half-array  inputs.  The  mean  and  variance  of  y.  the  correlator  out¬ 
put.  are  evaluated  as  follows: 

The  correlator  output  y(t)  is  the  product  of  the  derivative  of  the  process 
*^(t-T(t))  and  a  delayed  version  of  the  process  Xgft). 

y<t)  -  xj(t  -  r (t) )  *2  (t  -  ?(t>) 

Since 

*j(t>  =  s(t)  ♦  n ^(t) 

and 

*2(t>  a  •< t)  ♦  n2(t) 

where 


s.  n j .  and  n,  are  independent  zero  mean  gaussian  random  processes, 
it  is  straightforward  to  show  that  the  mean  of  y(t)  is  given  by 
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where 

o(t)  » 


t(t) 

t(t> 


normalized  autocorrelation  function  of  the  signal  or  the 
noise  processes 

time  value  of  the  moving  delay  parameter 
estimated  value  of  the  delay 
signal,  noise  power. 
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OUTPUT 


The  first  order  tracking  logic  updates  ?(t)  linearly  with  the  value  of 
y( t) .  On  the  average  when  t  ( t)  =r(t).  then  o’  lTlt)  -  t(t))  *  0  and  no 
changes  are  made  to  the  delay  estimate.  Thus 

r(n  ♦  1)  =  tin)  ♦  a  yin) 

The  process  yin)  can  be  represented  as  the  mean  value,  E  [y(n)]  ,  given 
in  (G-8)  plus  a  zero  mean  fluctuating  term.  N(n).  so  that 

tin  ♦  1)  =  T(n)  ♦  a  [E  (yin))  ♦  N(n)J 

2  *  A 

=  ?(n)  ♦  a  (-  !Jg“  o  lx(n)  -T(n))  ♦  N(n)l 
The  error  in  the  delay  estimates.  e(n).  is  given  by 
c(n)  «  T(n)  -  r (n) 

so  that 

“(n  ♦  1)  * *(n)  ♦  a(-  og2  o'  (f(n))  ♦  N(n)l 

Expanding  o  (c(n))  about  the  tracking  solution,  i.e.,  c(n)  *  0.  and 
using  only  the  first  term  in  the  resulting  Taylor  series, 

1  ’  2  '*  . 

-  <j  o  (c(n))  a  -  a  c(n)  o  (o) 

s  “ 

and 

r(n  ♦  1)  *  r(n)  ♦  a[-cg2  c(n)  o  (o)  ♦  N(n)J 

This  suggests  the  linearized  mathematical  model  for  the  split  beam  correlator 
tracker  shown  in  Figure  G- 2.  The  tracking  logic  is  modelled  as  an  integrator 
and  a  gain. 

G-8 


N(«) 


The  power  of  the  noise  process  NCn)  is  given  by  the  variance  of  y(n) 
and  is  assumed  to  be  independent  from  sample-to-sample.  Thus  it  remains 
to  compute  Var  l y ( n ) )  to  fully  specify  the  first  order  linearized  model. 

E  (y'(n)l  =  E  •{  l  s '( t-x)  ♦  n^’U-T)]*  I  s(  t-T)  *  n2(  t- ?)  1 2  J1 


=  E  [(s'(t-T))2  s2(t-?)l  -  2ag2on2o”(o)  -onV(o) 


I  It  A  »  J  1  »»  A  *» 

=  -<Jo  (o)  ♦  2a  [o  (t-?)1  *  -  2o  o  “  o  (o)  -  i  o  (o) 
s  s  s  n  n 


The  variance  of  y(n)  is  found  by  subtracting  the  square  of  (G-8). 

2  2 

Var  ly(n)J  *  -  (a  *  ♦  a  2)  o  (o)  ♦  e  “  [p  (r-?)J 

s  n  s 

For  small  errors,  i.e.  when  tracking,  the  second  term  is  small  and 

2 

Var  (y(n)J  =  Var  (N(n)J  *  -  (a  “  ♦  o  2)  o  (o)  (G-9) 

s  n 

Since  the  system  in  Figure  G-2  Is  linear,  it  is  straightforward  to  compute 
the  errors  due  to  noise  and  the  errors  due  to  the  dynamics  of  x(n). 

The  z-transfer  function  from  i(n)  to  c(n)  is  given  by 


H(.»  .  £$  ’ 


1  -  Z 


1  ♦  (oK-1)  Z 


(G- 10) 


G- 10 


The  8-transfer  function  from  N(n)  to  c(n)  is  given  by 


G(s) 


S 


3  2 


l  ♦  UK-1) 


Thus,  the  fluctuation  error  variance  is  given  by  13] 


(G- 11) 


[c2(n)] 


=  >  VjrS(n)  tl  Gd) 


Git'1) 


-1 


dz 


a2  Var  N(n) 
“  aRT2  -  aKT 


(G- 12) 


where  Eq  t G - 12)  has  been  obtained  by  the  calculus  of  residues. 


Dynamical  Response  to  Ramp  Inputs 

In  order  to  make  a  meaningful  comparison  to  the  adaptive  filter  bearing 
tracker,  it  is  necessary  to  select  the  same  inputs  to  both  devices.  The 
adaptive  filter  has  been  studied  for  a  linearly  varying  time  delay.  Hence. 
T(nt)  =  cnT  where  T  =  sample  time,  and  the  z-transform  of  the  Input  is 


(G- 13) 


and 


c  (i)  *  t(z)HO)  * 


_ cT _ 

(s  -  (  1  -  aK ) )  [t  -  ll 


(G-14) 
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G-12 


It  is  now  assumed  that  the  denominator  of  Eq  (G-17)  is  given  by  the  slope  of 
the  mean  weights  in  the  neighborhood  of  the  peak  (assuming  the  value  of  the 
slope  for  the  lag  given  by  Eq  (G-4)  does  not  differ  significantly  from  the 
value  of  the  slope  at  the  true  value  of  delay).  It  is  also  assumed  that  the 
shape  of  the  envelope  of  the  weights  is  well  approximated  by  a  decay  of  the 
form  o5. 

Then 


JL  -nhcn 

3*  H  )  t  J 


1  -  <  1  -  u a  2)oc6 
n 


(i  -0«12 


Using  Eqs  ( G- 18)  and  (G-19)  in  Eq  (G-17)  yields 


(G- 19) 


G-14 


Note,  m  a  check,  that  when  c  -  0  and  o  ■  e 


— 


b 


Var*  t 


which  agrees  with  (1). 


Note  also  that  Eq  (G-20)  has  the  correct  behavior  with 


2 

uon  (if  c  is  small)(and  o. 

2 

If  c  is  sufficiently  large,  then  decreasing  jon  degrades  the  performance 
because  the  weight  adjustments  cannot  effectively  keep  up  with  the  input. 

That  is.  there  is  lag  (Eq  (G-4))  in  the  location  of  the  peak  value  of  the 
mean  weights  and  the  amplitude  of  the  mean  weights  also  decreases 
(Eq  G-19).  Ignoring  the  lag.  the  variance  of  the  estimate  of  the  peak 
increases  due  to  the  amplitude  effect. 

The  bearing  tracker  performance  can  be  computed  directly  from 
Eq  (G-12)  since  Var  c(n)  is  Var  x(n).  Thus,  using  Eqs  (G-9)  and  (G- 12)  and 
Figure  G-Q, 

-IT—,  ( G- 21) 

o  (o)l 


G-15 


Equating  the  dynamical  errors  of  the  two  systems  as  given  by  Eq  (G-5) 
(multiplied  by  6  to  yield  time  delay)  and  Eq  (G-16)  (with  T  *  «  )  yields 


aK 


bed 


Tnnr 


(G-22) 


By  using  Eq  (G-22)  in  Eq  (G-21),  assuming  bcd«  and  using  Eq  (G-20),  the 
ratio  of  standard  deviations  is  given  by 


Var*  t 


Var^r 


adaptive  filter 


o.6y 


bearing  tracker 


oc6  ♦  k£ 


,C6 


uo 


(G-23) 


which  can  be  Airther  approximated  as 


(G-24) 


Conclusions 

The  Adaptive  Filter  Bearing  Tracker  and  the  Split  Beam  Bearing  Tracker 
are  compared  on  the  basis  of  tracking  error  variance  for  the  same  dynamical 
bias  error  when  tracking  the  linearly  time-varying  delay  of  a  broadband 
signal,  in  steady-state,  both  trackers  display  bias  errors  that  are  propor- 
tional  to  the  rate  of  change  of  delay  and  Inversely  proportional  to  the  loop 
gains.  Expressions  are  presented  for  variances  of  each  device  which  dis¬ 
play  the  tradeoffs  between  slgnal-to-noise  ratio,  signal  bandwidth  amt 


G-16 


I 


linear  rate  of  change  of  delay.  The  relative  performance  of  the  two  devices 

is  shown  to  be  essentially  independent  of  input  signal- to- noise  ratio  and 

to  depend  primarily  on  the  ratio  of  the  product  of  the  tap  time  delay,  6 ,  the 

signal  bandwidth,  b.  and  the  rate  of  change  of  delay,  c.  divided  by  the 

2  2 

adaptive  filter  algorithm  gain.  ucn  .  When  bci«uon  «1.  the  adaptive 
filter  variance  is  0.69  times  as  large  as  that  of  the  split  beam  tracker. 
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ADDENDUM  Gt 

EVALUATION  OF  EQUATION  (G-l)  FOR  (l-uo  2)  <  pC<S 
- a _ 

Beginning  with  Eq.  (G-l)  and  splitting  the  sum  in  two  parts 


:k(n)] 


k£j/C 

W082  £ 
k=0 

n-1 

*u°e2  £ 

k>i- 


n-k-1 


(l-uon‘) 


(j-ck)6 


n-k-1 

2.  (ck-J)« 

(1  -  uon  )  p 


=  va  2  »<« 


k<J/c 


*  )  P 


c5 


n-k-1 


♦ 


,  .  2  c5 (n-1) 

U  <3  p 


n-1 

E 


k>j/c 


n-k-1 


The  first  summation  represents  a  convergent  power  series  since 

2  c4 

1  ■  uon  <  ^  an<*  0  <  !•  ^he  »«cond  summation  will  converge  if 


l-uo. 


n 


cS 


<  l 


Now  let 

l  =  I 


w  *  integer  such  that-^  >  l  but  -jj-  <  l  ♦  1 


(GI-1) 


G-l  8 


Then,  in  the  first  summation 


For  large  t  .  the  second  term  in  the  first  set  of  brackets  is  negligible  in 
comparison  to  one.  This  approximation  yields  the  expression  in  Eq.  (G2)  for 


}  *c  (n-1). 


For  jtc(n-l).  Eq.  (G-l)  reduces  to  a  single  summation 


•  a  a 

E[«j<n)]  »u<ts2  £  (l  -n<rn2> 


n-k-1 


,(|-ck> 


23!Vc«n.l)  V  (1-.._2>»C‘ 


UO  3  3 
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n-k-l 


2  ® 
a  U<3  1  0 
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[|-c<n-»]  '-['‘-'.'"“l 


2  c  & 

1  -  (1  -  »cn  )0 


(GI-5) 


Hence 


4  fj-c(n-l)] 


As  a  check,  note  that  at  J  *  c(n-i).  Eq.  (G3)  is  self  consistent. 

The  approximate  location  of  the  peak  weight  in  (GI-4)  and  (Gl-6)  can  be 
determined  by  using  the  approximation  J  *  cl  and  assuming  that  n-l-1  can  be 
approximated  by  a  continuous  variable. 


Then,  assuming  the  peak  occurs  for  J  •  c(n-l),  Eq.  (GI-4)  simplifier  to 


n-t-l 


(l-ucn‘) 


2  c  6 

1  -  (1  -uon  )p 
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.-cut 
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where 
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-8 


(n-l-i) 
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*  0 


c  S  -be  5 
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=  1  -  (l-ucn2) 

*  l  -  (l-uon2)  °"C< 


(Gl-9) 


Differentiating  Eq.  (GI-8)  on  x  and  setting  the  result  to  zero  to  find  the 
maximum,  yields 
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ADDENDUM  UII: 

EVALUATION  OF  EQ.  (Gl)  FOR  (1  -  uo  2)>pc<* 


Consider  first  the  case  of  c  (n-1). 


ksj/c 

- 

[WH 

2  JS  -cS(n-l)  V' 

•  ““s  a  0  2^ 

k=J 
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n-1  . 

2  TT*  2  n_K‘l  (ck-j)  s 

♦  U0#  /  a  (1  ~uo n  )  0  (GII-I) 

k  ?  J/c 

The  first  term  sums  in  the  same  manner  as  Eq.  (GI-2).  However,  because 

2  c* 

of  the  condition  (i-  >  o  .  the  second  term  sums  differently  than  in 
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Adaptive  Filters,  Adaptive  Bearing  Trackers,  Bearing  Accuracy.  Broadband 
Bearing  Trackers,  L.MS  Adaptive  Filters,  Linearly  Varying  Bearing  Narrow¬ 
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A  specific  split  beam  adaptive  trackers  using  the  LMS  adaptive  filter  is 
Investigated.  The  adaptive  tracker  does  not  require  a  priori  information 
on  input  statistics,  and  can  track  using  narrowband  and  broadband  energy 
simultaneously. 

The  narrowband  tracker  is  shown  to  perform  within  0. 5dB  of  the  Carmer- 
Rao  lower  bound  for  narrowband  targets.  The  bias  and  standard  deviation 


r;  ♦»  T»TT*  ■  o 


OD  Form  1473:  Report  Documentation  Page  (Continued) 

Block  20:  Abstract  (Continued) 

with  plane  wave  interference  present.  When  tracking  a  narrowband  input  the  biaa 
effects  due  to  Interference  dominate  in  areas  in  flucuations  in  most  cases.  It  is 
shown  that  the  adaptive  tracker  is  less  sensitive  to  changes  in  input  spectra  re¬ 
lative  to  those  assumed  a  priori  than  conventional  split  beam  trackers,  and  the 
performance  of  the  adaptive  tracker  for  targets  with  linerly  changing  bearing  is 
analyzed  and  compared  to  a  conventional  split  beam  tracker. 
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